UNCLASSIFIED 


ad  402  114 


(lepn&luced 
luf, the 

DEFENSE  DOCUMENTATION  CENTER 


FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION,  ALEXANDRIA,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 

Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


"AFCRL  -63-83 


r=i 

fH 

03 

O 


ASTIA  DOCUMENT  NO:  AD _ 

CONTRACT  AF  61(052)  -  536 


Technical  Repon 


ELECTROMAGNETIC  CONTROL  OF  THE  ATTITUDE  OF  SPACE  VEHICLES 


"0 


fvj  ~>2 

C  •“> 

O  ;'3 
(•*»**- 1  " 

{**• 

C/2 
*tf c 


Authors:  P.  Bruscaglioni,  A.  Consortini 
Project  Director:  G.  Toraldo  di  Francia 


CENTRO  DI  STUDIO  PER  LA  FISICA  DELLE  MICROONDE 
CONSIGLIO  NAZIONALE  DELLE  RICERCHE 


FIRENZE  -  ITALIA 


December  1962 


!T  * 


A  P  T  I  A 

.  t nrr  jnrf 


» b  *  ^  m  V  i 


The  research  reported  in  Fh is  document  has  been  sponsored  by 
Electronics  Research  Directorate,  AFCRL  United 'States  Air  Force, 
monitored  by  the  European  Office,  Office  of  Aerospace  Research. 


ASTI A  DOCUMENT  NO:  AD 


CONTRACT  A?  61(052)  -536 


Technical  Report 


ELECTROMAGNETIC  CONTROL  OP  THE  ATTITUDE  OF  SPACE  VEHICLES 


Authors:  P.  Brusca^lioni ,  A.  Consortia! 
Project  Director:  G.  Toraldo  di  Franc! a 


CENTRO  DI  STUDIO  PER  LA  FISICA  DELLE  MICROONDE 
CONSIGLIO  NAZI ON ALE  DELLE  RICERCHE 
FIRENZE  -  ITALIA 


December  1962 


The  research  reported  in  this  document  has  been  sponsored  by 
Electronics  Research  Directorate,  AFCRL  United  States  Air  Fore 
monitored  by  uhe  European  Office,  Office  of  Aerospace  Research 


SUMMARY 


A  theoretical  investigation  has  been  made  of 
the  possibility  to  apply  the  properties  of  elliptically 
polarized  waves  for  the  attitude  control  of  space  vehi¬ 
cles. 

Three  different  ways  of  applying  electroma¬ 
gnetic  radiation  has  been  discussed  « 

The  first  one  employes  an  antenna  placed  on 
a  space  vehicle  ,  and  radiating  with  circular  polariza¬ 
tion  .  A  reaction  torque  turns  out  to  be  applied  to  the 
vehicle « 

The  second  one  makes  use  of  the  the  torque 
experienced  by  a  receiving  antenna  ,  placed  on  the 
vehicle  ,  when  a  ground  station  radiates  towards  the 
vehicle  .  The  third  one  is  based  on  the  simultaneous  and 
tuned  emission  of  radiation  from  an  antenna  located  on 
the  vehicle  and  a  ground  antenna? 

A  comparison  between  the  torque  which  can  be 
obtained  and  the  torques  necessary  to  the  control  of  the 
attitude  of  a  space  vehicle  in  actual  situations  shows 
that  the  first  procedure  can  have  pratical  importance, 
while  the  torques  obtained  by  means  of  the  second  and 
third  procedures  seerr  to  be  too  weak. 
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INIBODIICTION 


The  control  of  the  attitude  of  a  space  vehicle  is  of 
great  importance  in.  many  cases  «  For  instance  high  directivity 
antennas  can  he  used  on  communications  satellites  only  if  the 
satellite,  or  the  parts  containing  the  antennas,  can  he  suita¬ 
bly  oriented. 

The  case  where  the  power  to  he  used  by  the  vehicle 
comes  from  solar  cells-  placed  on  paddles  to  be  directed 
towards  the  sun,  or  from  solar  thermal  devices,  which,  require 
large  mirrors  to  concentrate  light  on  an  absorber,  is  another 
example  of  the  necessity  of  attitude  control.  In  these  cases 
the  requested  accuracy  can  he  of  the  order  of  one  degree  . 
However  in  the  case  of  some  satellites  for  scientific  resear¬ 
ches  the  accuracy  required  is  much,  higher  .  For  example  the 
O.A.O.  project  provides  that  the  satellite  must  keep  its 
attitude  fixed  in  an  inertial  frame  of  reference  ,  within  an 
angle  of  a  tenth  of  second  of  arcs 

Two  kinds  of  problems  arise  :  the  problem  of  using 
a  suitable  sensor  in  order  to  get  information  about  the  at¬ 
titude  of  the  vehicle  ,  and  the  problem,  of  applying  a  torque 
to  the  vehicle  ,  in  order  to  correct  its  attitude. 

There  are  a  number  of  different  attitude  sensors,, 
based  on  the  use  of  photocells  sensitive  to  visible  or 
infrared  light  ,  of  gyroscopes,  of  magnetometers,  of  tele¬ 
scopes  capable  of  determining  the  direction  of  a  star  with 
a  great  accuracy  ,  On  can.  always  find  a  device  sufficiently 
accurate  for  his  purpose  . 


The  attitude  control  can  be  effected  by  several 
methods,  which  are  generally  divided  into  two  classes , namely 
by  passive  and  active  methods  >  In  the  first  case  the  control 
is  obtained  without  an  expense  of  power  ♦  An  example  of  this 
case  are  the  spin  stabilized  satellites.  A  passive  stabiliza¬ 
tion  can  be  also  obtained  by  means  of  the  actions  of  the  gra¬ 
vitational  gradient  on  vehicles  of  elongated  shape,  or  by  means 
of  the  earth's  magnetic  field  acting  on  magnetic  bars  located 
in  the  vehicle  . 

The  active  methods  consist  in  the  application  of 
a  torque  to  the  vehicle  with  expense  of  some  power. 

Many  methods  for  applying  a  torque  to  a  space  vehi- 
cle  have  been  used  or  proposed  such  as  (  )  :  emission  of  gas 
jets  ,  controlled  motion  of  indipendent  parts  of  the  vehicle 
(rotating  wheels),  plasma  emissions,  excitation  of  currents  in 
coils  interacting  with  the  earth* s  magnetic  field  ,  changes 
of  the  mass  distribution  within  the  vehicle  (effect  of  the 
gravitational  gradient) ,  use  of  movable  solar  sails  and 
rudders. 

The  torques  achievable  with  the  different  methods 
mentioned  turn  out  of  different  order  of  magnitude.  Thus  when 
a  quick  maneuver  of  a  vehicle  is  hecessary  ,  the  emission  of 


(^)  Bee  for  instance  Ref.1 


gas  jets  is  considered  to  be  the  most  suitable  device  because 
of  the  large  value  of  the  torque  obtainable.  In  other  cases 
however  very  weak  torques  are  sufficient;,  for  instance  when 
the  attitude  of  a  vehicle  must  be  kept  stable  against  the 
action  of  small  external  forces. 

The  object  of  this  research  is  a  discussion  of  a 
method  based  on  the  use  of  elliptic ally  polarized  electroma¬ 
gnetic  waves.  Such  waves  have  the  property  to  carry  angular 
momentum  about  the  direction  of  propagation. 

This  property  was  first  derived  from  theoretical 

considerations  by  Sadowsky  (  ),  Epstein  (“),  and  Poynting  (  ). 

5 

Later  it  received  experimental  proof  by  Beth  (  ),  who  used 

.  6. 

visible  light  ,  and  by  Carrara  (  )-  who  used  microwaves. 

Three  different  ways  of  employing  electromagnetic 
radiation  for  the  attitude  control  of  a  space  vehicle  will  be 
discussed  in  this  report,  and  will  form  the  subjects  of. the 
first,  second  and  third  sections  respectively.  The  first 
is  an  active  method,  the  second  a  passive  ro.ethod,  the  third 
a  mixed  method. 
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1  -  REACTION  TORQUE  ON  A  RADIATING-  ANTENNA 


1.1  -  General  remarks 

In  the  present  part  of  this  report  we  will  discuss  the 
possibility  of  employing  the  torque, experienced  by  an  antenna 
placed  on  a  space  vehicle, to  control  the  attitude  of  the  vehicle. 

First, a  formula  will  be  derived  for  the  calculation  of 
the  flux  of  angular  momentum  through  a  spherical  surface -This 
formula  will  then  be  applied  to  the  evaluation  of  the  torques 
acting  on  the  antennas.  In  the  sequel  such  torques  will  be  indi- 
cated  as  ' 1  electromagnetic  torques'1  . 

Consider  a  radiating  antenna, which  for  simplicity  will 
be  assumed  to  be  a  point  source, located  at  P.  It  is  possible  for 
the  electromagnetic  radiation  issued  by  the  antenna  to  have  an 
angular  momentum  about  P.  This  can  happen  when  the  radiation  is 
elliptlcally  polarized-*  In  this  case  it  is  evident  that  the  an¬ 
tenna  must  be  subjected  to  a  reaction  torque. 

1.2-  Flux  of  angular  momentum  through  a  spherical  surface 
containing  charges  and  currents  . 

Let  V  be  a  region  of  space  in  which  charges  and  cur¬ 
rents  are  present, and  £dV  the  momentum  of  the  charges  contained 
in  a  small  volume  dV.  Let  S  be  the  Poynting  vector  and  T  the 
Maxwell's  tensor  defined  as  : 

T-ED  +  HB  -  UU 
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where  E,  H  are  the  electric  and  magnetic  field  vectors,  D,  B 
are  the  electric  and  magnetic  induction  vectors.  IXis  the 
density  of  electromagnetic  energy  : 

U  =  \  (EcD  +  H,B) 


and  U  is  'the  fundamental  tensor  of  the  Euclidean  space, 
whose  components  are  : 


U 


lm 


-s 

"  lm 


L  =  0  if  1  /  i 
lm 


1  if  1  =  m 


By  taking  into  account  the  action  of  the  field  on  the  charges 
and  currents,  it  is  possible  to  obtain  the  following  rela¬ 
tion  (  1 )  : 


(1.2-1) 


9Q 

at 


+  div  T 


where  v  is  the  velocity  of  the  electromagnetic  waves  in  the 
medium.  Eor  the  angular  momentum  of  the  charges  with  respect 
to  a  certain  fixed  point  0,  one  obtains  from  (1,2-1): 


(1.2-2)  L  (QAr) 


a 

at 


(div  T) A  r 


where  r  denotes  the  radius  vector  of  the  point  of  interest  with 
respect  to  0. 


( ' )  see  for  instance :Ref„7  p.?8 
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The  total  angular  momentum  of  the  charges  will  he 


given  hy: 


(1,2-3) 


ht 


Qa£  dV  *  -  tj-  /  dV  +  j  (div  T)at  dY 

vy  \/y 


Since  the  following  relation  is  valid  for  a  dyadic  u  v 
(see  Appendix  1) 


(1  2-4) 


div  (u  v) 


A ' 


div 


(u_Ar) 


uav 


the  last  integral  in  (1,2-3)  can  be  written  as  : 


(1,2-5)  /  (div  T)at  dV  =  I  div 


V 


(EAr)  D' 


dV  + 


+  div 

h 


(Hat)  B  dY  +  I  (divUU)Ar  dV 


If  the  region  Y  is  spherical  with  its  center  at  0,  it;  is 
possible  to  demonstrate  that  the  last  term  of  Eq_.  (1.2-5) 
is  zero  (see  Appendix  2) »  Thus,  by  applying  the  divergence 
theorem,  one  has  : 


6 


(1.2-6) 


a 

“  3t 


f 

L'a.rav  . 

/  V 


SAr  dV 
v 


+  I  (E.n)  DAr  dl  +  (H.n)  Bat  <jE 

4  4 


where  n  is  the  outward  normal  to  the  surface  2  containing  V. 
The  expression: 


m  -  *• 


Sat  dV 


is  the  resulting  moment,  with  respect  to  0,  of  the  electroma¬ 
gnetic  forces  applied  to  the  charges  in  V. 

If  £>  H,  D,  ®  are  oscillating  functions  of  constant 
amplitude,  the  time  average  M  of  m  is: 


(1.2-7) 


M  = 


1 

2 


He  j  (B*.n)  DAr  dl  +  i  Re  I  (^n)  BAr  d2 


where  the  asterisk  x  denotes  the  complex  conjugate  quantity. 
The  time  average  of 

f 
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is  in.  fact  zero  .  Some  situations  c^n  occur  ,  in  which.  M  is 
different  from  zero  . 

1 . 3-Torque  on  two  crossed  dipoles  in  quadrature 

As  a  first  simple  example  let  us  consider  the  case 
in  which  the  antenna  is  formed  by  two  elementary  dipoxes  , 
perpendicular  to  each  other  (Pig.  1.-1 )  .The  currents,  assumed 
to  be  uniform  across  both  dipoles,  have  equal  amplitudes:  arid 

71 

their  phases  differ  by  ^  . 

The  antenna  is  assumed  to  be  placed  in  the  vacuum, 
far  from  other  bodies  .  Thus  the  field  of  the  antenna,  in  the 
region  considered  for  the  following  calculations,  can  be  assu¬ 
med  as  being  undisturbed  by  the  presence  of  other  bodies  • 


All  elements  of  dipole  1  will  be  labelled  by  the 
index  1,  all  elements  of  dipole  2,  by  the  index  2.  Dipole  1 


is  assumed  to  "be  oriented  along  the  z  axis  (unit  vector  Jf.)  of 
a  rectangular  coordinates  system,  -whose  center  is  placed  at 
the  point  0  where  the  dipoles  cross  each  other  «  Dipole  2 
is  oriented  along  the  y  axis  (unit  vector  * 

The  currents  of  the  two  dipoles  will  "be  represented 
"by 


I 


1 


*»  le 


i'f't 


(1.3-1) 


.  T  iwt 
X  le 


where  i  is  the  imaginary  unit  « 

let  us  consider  two  systems  of  polar  coordinates? 

T-} 5  and  r ,  Og,  cpg,  both  with  their  origin  at  0  and  with 
the  polar  axes  along  the  z  and  y  directions  respectively  (see 
Pig .1-2) -Let  £  represent  a  spherical  surface  with  center  at  0 
and  radius  r  »  The  antenna  radiates  with  circular  polarization 
in  the  direction  x  ,  The  reaction  torque  applied  to  the  antenna 
can  be  evaluated,  by  means  of  (1.2-7)  »  by  calculating  the  flux: 
of  angular  momentum  through  the  surface  £  * 

It  is  seen  at  once  that  in  this  case  the  expression  : 


1 

"2 


Re  ji 


(Hs.n)HAr  d£ 

4, 


which  appears  in  (1,2-7)  vanishes,  because®  the  magnetic 
field  vector  is  tangent  to  the  surface  £  «  Thus  (1,2-7)  can  he 
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Polar  coordinates  r,  ©1  ,  cp1  and  r,  0 2>  cp2> 

n,  i  ,  i  ,  are  the  unit  vectors  of  the  reference  system  1,  and 
Q1  cp  1 

n,  ir.n»  i_o  are  the  unit  vectors  of  the  reference  system  2, 


written  as  : 


(1.3-2) 


He  f. 


o 


(ifhnj  2Ar  dl 


k 


and  the  component  of  M  along  the  x  axis  is; 


(i«3~3) 


M 

x 


-  Re  e 
2  o 


(E  ,n)  E/.r.i  d£ 


If  E  and  E  represent  the  transversal  and  radial  components 

'ts  T 

of  E,  and  s  denotes  the  common  length  of  the  dipoles,  one  has, 
apart  from  terms  of  the  order  l/r^, 


-1  ”1r  ~  +  E10  "01 


(1.3-4) 


E 


1r 


E 


Z  Is 

c  S  cos0i  e 

r 


10 


.  Z  Is  n 
x  -  T-  senO„  e 
2  hv  1 


i(wt-kr) 

i(mt-hr) 


-  11  - 


—2 

E2r  a  +  E20  ±Q2 

(1.3-5) 

E2r 

iZ  Is  n  i(wt-kr) 

‘  K  ~l  003°£  e 

r 

E20 

Z  Is  A  i(tot-kr) 

’-■5  1?  sanS2  e 

Putting  E 

-  Si  • 

f  E2  ,  one  has: 

(1.3-6) 

E*.n 

-  (®i  +  52>-e  -  E,r  +  E2r 

besides  : 

(1.3-7) 

Ea  n 

*  (I,  +  S2'as  =  *  Eia  i^i  ■  E2a  3^2  > 

X 

=  r  n 

Thus  (1.3-3)  becomes: 

(1.3-8) 

K  = 

X 

-  r  Ee  j(E*.  +  E2r><E1«V  +  E2«i<p2)'ir  41 

From  Fig .1-2  one  has: 
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(1.3-9) 


i^.i  =  -  sincp. 


i^.x  =  -  sxncp2 


Hence,  with  the  substitution  &E  =  r  dQ,  (1.3-8)  becomes: 


e  r 
o 


(1.3-10)  R el  (E*r  +  E^r)(Eiesincp1  +  E^sintpg) 


dQ 


where  the  integration  is  over  the  whole  solid  angle  4n. 
By  means  of  Eq..(  1.3-4)  and  (1.3-5)  one  has: 


e  Z 


2  2/ 

(1.3p11)  =  ^TiX  *  3  (-cosO^sinOgSincpg  +  sinO^ cos^sincp  1 


)dQ 


Prom  Pig. 1-2  one  cam  find 


coatf^  =  sirup^sinO^ 

(1.3-12) 

cos01  -  -  sin(p2sin02 


-  13  - 


Hence  ( 1  *  3— 11)  "becomes  : 


e  Z 


1 2n 


'% 


(1.3-13)  l"s2  I  sin2^  dtp  !  sin3«  d .Q 

'o  1 0 


x 


4  £0Z  2  2  1 

3lTis  =5  3  % 


t2  2 
-i-  B 


It  is  seen  that  the  components  and  M  of  M  are  equal  to  zero, 
since  their  expressions,  corresponding  to  the  ( 1  .,3-1 3),  are  of 
the  form: 


'2* 


r 


<„3, 


sincp  coscp  dcp  i  sin-7  3  d$ 


/° 


The  power  W  radiated  "by  the  two  dipoles  is  gi9en  by: 
(1.3-14-)  W  =  W1  +  W2  =:  1  u 


By  comparison  with  (1.3-13)  o:ae  obtains  the  relation  between  the 
radiated  power  and  the  torcue  applied  by  reaction  to  the 
antenna: 


( 1 .3-15) 


X  •-= 


W 

5 


The  torque  M  is  directed  along  the  normal  to  the  plane  of  the 


14  - 


dipole 3. > 


T'he  antenna  tends  to  rotate  in  a  direction  opposite 
to  that  of  the  positive  (or  negative)  charges  at  the  end  of 
the  dipole':. 

It  is  interesting  to  note  that  result  (1.3-15)  is 
exactly  the  same  as  one  would  get  in  the  case  of  a  plane  wave 
with  circular  polarization  ,  although  in  our  case  the  polari¬ 
zation  is  perfectly  circular  only  in  the  x  direction  .  The 
apparent  paradox  can  he  explained  by  considering  that  in  our 
case  a  non  vanishing  contribution  to  angular  momentum  comes 
from  the  radial  component  of  E,  which  combined  with  the 
tangential,  component  of  H,  gives  rise  to  a  tangential  momen¬ 
tum  . 


1*4-  Torque  acting  on  a.  turnstile  antenna. 

The  calculation  of  the  preceding  section  will  now  be 
extended  to  the  case  of  a  turnstile  antenna  placed  in  the  va*- 
cuum  sind  far  from  other  material  bodies  .  The  antenna  is 
constituted  by  two  half-wave  dipoles  fed  with,  currents  of  equal 
amplitudes  and  in  quadrature  .  The  eirms  of  the  eintenna  are  as¬ 
sumed  to  be  infinitely  thin,  and  the  currents  to  have  a  sinusoi¬ 
dal  distribution  .  With  reference  to  Pig*3«>3  we  will  write: 
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Pig. 1-3 

Reference  systems  x,  y,  z;  r,  0^  cp  ^  *  r,  ©2>  cp2  for  the 
turnstile  antenna,  formed  by  the  two  half-wave  antennas 


1  and  2. 


(1.4-1) 


Ij  »  X  oos(kz)eiu,t 

Ig  *  i  I  cos(ky)eiu1:  ,  k |- 


with  a  notation  analoguos  to  that  of  the  preceding  sections. 
Eq.  (1.2-7)  will  be  used  for  the  evaluation  of  the  torque 
acting  on  the  antenna. 

The  expressions  of  the  electric  fields  are  (*) 


Si  - 

E1l^  *  Eici01 

(1.4-2) 

E1r  “ 

---«  3in(^  cosO 
8*r  V2  1 

f  \ 

j  ei(wt-kr) 

1ZI  c°allcoaai) 

i(ut-kr) 

EiB* 

2.5  SI557 

e 

E2  “  E23~  +  E2fii02 


(1.4-3) 


,  s  -21 
J20 


y 

oob(|oo3oJ  i(ut.ltr) 


_  iZXX  .  [n  „„  n  \  i(ut-kr) 

hr  *  ~ 5  alnb  003°W* 


’sln0~ 


(')  See  Ref .8,  p.528 
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With: 


E  =  E.  +  E_ 
—  —1  —2 


E*.n  =  E*  +  E* 
—  —  1r  2r 


Eaii  =  -  i  „E  ~  i  0Ert/N 

-  — rp  1  10  — cp  2  2  0 


Eq_.  (1,2-7)  becomes: 


0.4-4)  S  -  -  5*  *•  /  <*fr  +  *  S24^2)r  « 


because  the  magnetic  field  vectors  are  tangential. 
Substitution  of  (1.4-2)  and  (1.4-3)  into  (1.4-4)  gives: 


or  also : 
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A 

1  sin(|  cos© 2]  cos! 

[l  c0,,il 

sinj 

„  ±  _ 

|  cosQ.jjcos 

(f  00B°^ 

ii53~' 

V 

iinfl" 

- 

z 

2 

By  taking  into  account  (1.3*9),  with  the  substitution  d£  =  r  dQ 
one  obtains: 


(1.4-7) 

L 


£  Z  „ 

M.i  »  I  X  x 

32n 


y 


sintp^in^cosOg^cos^008®^  sirupgSin  ^cosoj  cosj?|cos9^ 


'slnS: 


1153; 


dQ 


where  the  integration  is  over  the  whole  solid  angle  Q  ■  4ij. 
By  introducing  (1.3-12),  Eq.(  1.4-7)  becomes: 


(1.4-8) 


e  Z  „ 

-2-,  I2  * 
16* 


sinq>  Bin(j  sincp  sinojcos  (joos©) 
iln3 


dQ 


dZ 


which  gives: 


( 1  • 4—9 )  It  -  — I  X  2ti 

x  1 6n 

i 

By  using  the  following  formula,  whose  demonstration  will  he 
given  in  Appendix  3  , 


(1.4-10) 


r 

*inO)  cos(P  cos Q)  dO  = 

lo 


from  (1.4—9)  one  has  (a  =  p  =  %/2): 


-  0.605699 

u 


Hence  e<i. (1.4—9)  becomes; 


(1.4-11) 


M 

x 


I2X  0.605699  =  g-  |  I2  0.605699 
2 n  C  n 


By  introducing  the  radiation  resistance  of.  a  half-wave  component, 
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given  "by  (*  )  * 


(1.4-12)  R  =  |j  2.43765 

and  the  power  radiated  by  the  turnstile  antenna,  one 

2 

obtains  recalling  that  w  ■  I  R, 


(1.4-13) 


X  =  0.9939  - 

x  w 


The  calculation  of  the  torque  applied  to  the  antenna 
has  been  performed  also  in  another  way  ,  namely 4  by  talcing 
into  account  the  mutual  actions  between  the  currents  and  the 
charges  on  the  four  arms  of  the  antenna.  This  second  way  is 
longer  and  more  complicated  than  the  one  described  above,  and 
the  calculations  are  shown  in  Appendix  4 .  Here  only  the  results 
are  given,  which  were  obtained  by  the  use  of  some  approximations. 

The  value  for  M  is 
x 

M  =  ‘  0.98  - 

x  to 


('•)  see  for  instance  Ref.  8  p.148  Eq.(l6) 


Though,  the  result  is  not  exact  ,  this  procedure  may 
he  more  convincing  since  it  is  a  direct  physical  approach  . 

It  is  seen  at  once,  hy  considering  the  mutual  forces 
between  charges  and  currents,  that  the  x  and  y  components  of  M 
are  equal  to  zero  . 

1.5  -  Application  to  the  attitude  control  of  space  vehicles. 

Eq*  (1.4-13)  of  the  preceding  section  shows  that  the 
reaction  torque  applied  to  a  turnstile  antenna ? assumed  to  he 
isolated  in  the  vacuum,  is  proportional  to  the  radiated  power 
and  inversely  proportional  to  the  frequency  .  If  a  turnstile  an¬ 
tenna  is  placed  on  a  vehicle,  as  shown  in  Pig^l-4.  the  whole 
vehicle  is  submitted  to  a  torque  ,  when  the  antenna  radiates. 


Ei g,1-4 
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Let  us  see  how  this  source  of  torque  lends  itself  to 
he  applied  to  the  control  of  the  attitude  of  the  vehicle* 

First  the  value  of  the  torque  which  can  be  obtained  ,  will  be 
considered  . 

Eq.  ( 1  3)  was  derived  with  the  assumption  that 

the  far  field  of  xne  turnstile  antenna  was  undirsturbed  by 
the  presence  of  other  bodies.  If  the  size  of  the  vehicle  is  so 
small, in  comparison  with  the  antenna,  that  the  far  field  of  the 
antenna  is  not  greatly  modified,  the  value  of  the  torque  can  be 
assumed  to  be  given  by  (1.4-13). 

In  order  to  obtain  numerical  values,  one  must  know 
the  level  of  the  electric  power  available  on  a  space  vehicle  * 

Up  to  now,  solar  cells  or  chemical  batteries  have 
been  used  as  power  sources  on  space  vehicles.  In  one  case  a 
nuclear  thermo-electric  generator  has  been  employed. 

From  the  data  reported  in  the  literature  ,  it  seems 
that  the  case  in  which  the  highest  level  of  power  was  used,  is 
that  of  the  Courier  (I960  v)  1  B  :  225  watts  were  available  for 
the  transmission  to  ground  stations.  The  power  was  supplied 
by  storage  batteries  ,  charged  by  a  system  of  19152  solar  cells, 
capable  of  a  continuous  power  of  60  watts. 

For  the  Nimbus,  the  meteorological  satellite  which 
Was  scheduled  for  the  end  of  1962,  the  power  available  from 
batteries*  charged  by  solar  cells,  was  planned  to  be  250  watts. 

In  the  other  cases  reported  in  the  literature,  at 
most  some  tens  of .  wktts’  where  used. 

The  forecast  is  that  ,  in  the. future  ,  the  electric 


power  necessary  for  space  missions  will  be  much  higher.  A 
communication  satellite  for  broadcast  service  would  need 


approximately  1500  watts  of  power  .  50  kW  of  power  have 

1 0 

been  indicated  for  a  television  transmitter  (  )  , 

11 

Snyder  (  )  indicates  that  tens  of  kW  will  be  neces¬ 

sary  for  Space  Stations  , hundreds  of  kW  for  manned  satellites 
of  long  life,  and  up  to  20  MW  in  some  cases  of  big  manned  and 
waaeorr  rable  v eb.tcl  e  s . 

Therefore  several  methods  of  producing  electric 
power  are  under  study,  in  addition  to  those  which  have  been 
already  applied. 

Solar  cells  could  give  power  up  to  some  kW,  though 
they  are  not  considered  convenient  for  power  higher  than  some 
hundreds  of  watts. 

Solar  thermoelectric  and  thermoionic  system  are  being 
studied,  where  a  mirror  concentrates  solar  energy  on  an  ab  - 
sorber,  in  order  to  heat  a  junction  of  a  thermocouple  ,  or  the 
cathode  of  a  thermoionic  valve. 

With  these  systems,  it  seem  possible  to  achieve  power 
levels  up  to  25  KW„ 

Solar  energy  can  be  also  employed  to  feed  a  turbo¬ 
engine  by  means  of  mirrors  concentrating  the  solar  heat.  This 
system  is  indicated  for  power  of  some  tens  of  kW,  As  ah  example 

of  the  weight  of  these  system^,  one  may  mention  a  pro ject , under 

1 2 

study  at  the  Sunstrand  Denver  Corporation  (  )  ,  for  a  solar 

turbo electric  engine  of  15  kW  ,  having  a  weight  of  820  lb  , 

It  is  forseen  that  such  a  weight/power  ratio  of  55  lb/fcW,  can 
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be  lowered  to  a  value  of  20  Uo/kW ,  by  improvement  of  the  tech¬ 
nique  . 

Nuclear  energy  is  considered  suitable  for  producing 
electrical  power  on  a  space  vehicle.  The  satellite  Transit  IV  A 
used  a  source  of  power ,  consisting  of  a  thermocouple,  whose 
hot  junction  was  placed  in  a  small  reactor  containing  Pt  238. 
This  device  delivered  2.7  watts  and  its  weight  was  2.1  Kg. 

Several  systems  based  on  nuclear  energy  are  under 
study.  The  availability  of  the  three  units:  SNAP  10,  SNAP  2, 
SNAP  8  is  expected  for  the  years  1963-65.  The  three  systems., 
will  provide  300  watts,  3000  watts, 35000  watts,  respectively, 
with  weights  of  the  units  of  300,  600  and  1500  lb  ( 1 . 

The  availability  of  the  system  SNAP  50,  delivering 

*1  A 

300  kW,  is  forseen  for  the  year  1970  (  )  . 

Prom  Eq. (1.4-13)  it  is  seen  that  the  torque,  which 

can  be  obtained  by  placing  a  turnstile  antenna  on  a  vehicle, 

is  proportional  to  the  wavelength  A  .  If  a  good  electrical 

efficiency  is  desired  ,  an  upper  limit  for  A  is  set  by  the 

dimensions  of  the  antennas  which  can  be  placed  on  the  vehicle. 

The  longest  antenna  placed  on  a  space  vehicle,  up  to  now,  is 

the  150  ft.  radio  antenna  of  the  Canadian  satellite  Alouet- 
15 

te  (  )  .  Generally  the  dimensions  of  the  satellite  an¬ 

tennas  were  much  smaller  ,  of  the  order  of  one  meter. 

Table  I  ,  calculated  by  means  of  (1.4-13)  ,  gives 
the  torques  which  can  be  obtained,  with  power  up  to  some  kW.? 
and  antennas  capable  of  radiating  with  wavelenghts  up  to 
100  m. 
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TABLE  I 


Values  of  the  electromagnetic  torques  E 


M  =  5.3  W  A  ,  where  Iff  is  given  in  ergs,W  in  kW,A  in  meters. 


W  (kW) 

A  -  5  m 

M  (ergs) 

X  =  10  in 

A  =  20m 

A  =  100  m 

0.1 

2  .65 

5-3 

10.6 

53 

1 

26  v  5 

53 

106 

530 

5 

1 

132.5 

265 

530 

2650 

The  torques  requested  for  the  attitude  control  of 
the  vehicle  vary  largely  according  to  the  type  of  vehicle, 
the  mission  to  be  performed  and  the  particular  period  of  its 
activity  •. 

First  the  principal  causes  ,  which  can  disturb  the 
attitude  of  a  stabilized  vehicle  ,  will  be  considered  . 

The  principal  disturbing  actions  to  be  taken  into 
account  are  the  following  : 

Gradient  of  the  gravity 
Inertial  reaction  forces 
Earth's  magnetic  field 
Solar  radiation  pressure 
Meteorites  impact 
Atmospheric  drag 
Motion  of  internal  parts 
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These  actions  will  be  discussed  one  by  one,  in  the 
following  sections. 


1.6-  Gravitational  gradient  and  inertial  reaction  torques 

a)  Gravitational  gradient. 

A  torque  is  exerted  on  a  vehicle  of  non  spherical 

mass- symmetry  in  the  proximity  of  bodies  of  great  mass.  This 

action  occurs  when  the  axis  of  minimum  moment  of  inertia  of  the 

vehicle  is  not  directed  towards  the  center  of  gravity  of  the 

attracting  body.  In  the  case  of  an  earth's  satellite,  this 

gives  rise  to  torques  G  and  G  about  the  pitch  and  roll  axes  ( r ) 

P  r1 5 

and  their  values  are  given  by  (  ) 


G  =  i  u)2  AK  sin 20 
p  2  o  p  p 


(1.6-1) 


G  =  ^  u2  AK  sin 20. 
r  2  o  r  r 


where  AK  , 
P 


AK  are  the  differences  between  the  moments  of 
r 


(,•')  For  earth  centered  references*  generally  pitch  axis  indicates 
the  axis  passing  through  the  center  of  mass  of  the  satellite  and 
normal  to  the  orbital  plane  .  The  yaw  axis  in  the  axis  directed 
along  the  vertical  direction. The  roll  axis  is  the  axis  normal 
to  the  preceding  two. 
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inertia  about  the  pitch  or  roll  axes  and  the  moment  of  inertia 

about  the  vertical  direction.  0  and  0  are  the  attitude  de- 

P  ^ 

viation  angles  from  the  stable  direction,  in  pitch  or  roll  . 

is  the  mean  angular  orbital  velocity  . 

In  the  case  9  =  0  =  45°  and  AK  =  AK  =  AK, 

p  r  p  r 

the  expression  G/AK  =  G  /AK  =  G  /AK  has  been  evaluated,  a3 

p  p  r  r 

a  function  of  the  altitude  .  The  results  are  reported  in 
Table  II . 


TABLE  II  -  Gravitational  torque  G/AK 


Altitude  ( km) 

G/AK  (j/kg  m2) 

200 

2.1  10“6 

400 

1.9  10“° 

600 

1.8  10~6 

1000 

1.5  10”6 

2000 

1  io“6 

3000 

0.64  10“6 

4000 

0.53  10"6 

5000 

0.36  10"6 

10000 

1.3  10'7 

15000 

0.6  10"7 

The  effect  of  the  gradient  of  the  gravity  on  an 
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elongated  vehicle  can  be  employed  aa  a  means  for  stabilizing 
the  attitude  (example:  the  Iraac  satellite)  .  In  general  it 
represents  a  cause  of  disturbance  for  a  stabilized  satellite, 
and  must  be  neutralized  .  Prom  Table  II,  it  is  seen  that  the 
torque  to  be  balanced  can  be  of  the  same  order  of  magnitude  as 
the  torques  obtained  by  means  of  emission  of  oircularly  polarized 
radiation. 

b)  Inertial  reaction  torques. 

In  the  studies  on  the  attitude  control  of  space 
vehicles,  the  consideration  of  the  particular  reference  sy¬ 
stem  is  of  primary  importance  ( ^ )  .  Por  example  an  inertial 
system  is  suitable  for  describing  the  attitude  of  space 
astrophysical  observatories,  while  a  geocentric  system  is  ge¬ 
nerally  used  for  communications  or  meteorological  satellites. 

Let  us  suppose  that  a  space  vehicle  must  keep  its 
orientation  fixed  with  respect  to  a  non-inertial  reference 
frame  .  Then  the  equations,  describing  the  motion  of  the 
vehicle  in  this  reference  system  ,  contain  some  terms,  which 
represent  torques  applied  to  the  vehicle  ,  due  to  inertial 
reactions. 

As  an  example  of  these  Inertial  reaotion  torques, 
let  us  consider  a  communications  satellite  ,  which  must  keep 
an  axis  (e.g.  the  maximum  directivity  axis  of  an  antenna)  di¬ 
rected  towards  the  earth  .  In  this  case  it  is  useful  to  con¬ 
sider  a  rotating  reference  frame  ,  with  angular  velocity  equal 
to  the  angular  orbital  velocity  of  the  satellite.  If  the  orbit 
is  elliptic  ,  with  eccentricity  e,  the  inertial  torque  T, 
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-1 

acting  on  the  satellite  is  approximately  given  by,(  ) 


(1.6-2)  T  =  2  K  e  y  sinu)  t  =  T  sinu)  t 

Poo  oo 

v/here  Kp  is  the  moment  of  inertia  of  the  satellite  ,  with 
respect  to  the  pitch  axis  ,  is  the  mean  orbital  angular 
velocity. 

The  direction  of  the  torque  is  normal  to  the  plane 
of  the  orbit. 

We  will  evaluate  T  in  some  practical  cases. 

1.7  -  Some  examples: the  synchronous  satellite,  satellites 
orbiting  at  other  altitudes 

Let  us  consider  a  synchronous  communications  satel¬ 
lite  ,  orbiting  with  a  period  of  24  hours  at  a  height  of 
~  36000  km.  Prom  this  height  the  earth  is  seen  under  an  angle 
of  v17°  .  It  is  therefore  advisable  to  use  a  good  directivity 
antenna  ,  directed  towards  the  earth. 

The  torque  T  is  equal  to  zero  ,  when  the  orbit  is 
circular  .  In  this  case  the  desired  attitude  is  maintained,  if 
the  satellite  has  been  placed  in  orbit  with  a  spin  velocity 

=  u)  about  the  pitch  axis  . 

0 

However  ,  when  the  orbit  is  elliptic,  a  torque  must 
be  applied  to  the  satellite  to  counterbalance  the  inertial 
reaction  torque  ^ ' 

Let  us  consider  a  parallelepiped  satellite  ,  with  the 


-  30  - 


geometrical  situation  described  in  Fig ,1-5.,  It  is  convenient  that 
the  sf>in  axis  is  the  axis  of  maximum  moment  of  inertia  and  that 
the  axis  of  minimum  moment  of  inertia  is  directed  towards 
the  earth,  so  that  the  attitude  is  inherently  stable  . 


i 


l 

i 


Fig. 1-5 
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Table  III  gives  the  maximum  value  T  of  the  inertial  tor¬ 
que  , calculated  by  means  of  (1*6-2).  Several  values  of  the  eccentri¬ 
city,  of  the  size  and  weight  of  the  satellite  are  considered. 


TABLE  III 


Maximum  torques  on  the  satellite  of  Pig. 5,  due  to  inertial  reaction 
force.  The  sides  a,  b,  c  are  given  in  m,  the  weight  in  kg,  the  tor¬ 
que  in  J.  The  eccentricity  is  denoted  by  e. 


i 

i 

1 

c 

(m) 

weight 

(kg) 

.  X 

P  2 
(leg  m  ) 

<D 

11 

°, 

e=10 

Torque  Tq  (J) 

-2  <Hr1 
6*10^ 

e*2.10-1 

1 

1 

2 

500 

208 

2.21 

io-9 

2.21 

10-6 

2.21  ,icf7 

1 

1 

2 

1000 

417 

4.42 

IO"9 

4.42 

■  10“8 

4.42-10”7 

8.84 

io-7 

1 

1 

2 

2000 

834 

CO 

« 

oo 

4* 

_q 

.10  * 

8.84 

10-8 

8.84  10“7 

1 .76 

io”6 

1 

1 

3 

1000 

834 

8.84 

10"9 

8.84 

,o-8 

8.84- 10“7 

1 .76 

io-6 

1 

1 

3 

2000 

1665 

1.76 

,o-8 

1.76 

10-7 

1.76  10~6 

3.52 

io-6 

2 

1 

5000 

8330 

00 

• 

oo 

VJ1 

io”8 

8.85 

10-7 

8.85  10**8 

1 .77 

io-5 

2 

B 

6500 

10800 

1.15 

10"7 

1.15 

10-6 

1.15  10“5 

2.3 

io-5 

2 

B 

8000 

13320 

1 .41 

10-7 

1.41 

io'6 

ir\ 

I 

O 

T* 

• 

2.8 

io-5 

2 

H 

10000 

16700 

1.77 

10-7 

1.77 

10-6 

1.77  10“5 

3.54 

10“  5 
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If  the  axis  of  minimum,  moment  of  inertia  makes  an 
angle  0  with  the  direction  of  the  gravity  ,  the  action  of  the 
gradient  of  the  gravity  results  in  a  torque  ,  applied  to  the 
satellite,  which  tends  to  lower  0  .  Thus  if  the  inertial 
torque  tends  to  make  the  satellite  deviate  from  its  attitude, 
a  restoring  torque  G  arises  from  the  action  of  the  gradient 
of  the  gravity  .  Precisely  G  is  given  by  : 


(1.7-1) 


G  =  -  }  W2  (K  -K  )sin2© 
2  o  p  y 


where  K  and  K  are  the  moments  of  inertia  with  respect  to  the 

p  y 

roll  and  yaw  axes  respectively  (Fig.1-5). 

If  other  disturbances  can  be  neglected  ,  and  there 
is  no  attitude  control,  the  equation  for  Q  can  be  written  as  : 


(1.7-2) 


K  0  »  2  K  e  w2  sinu  t  -  |  u2  (K  -  K  )sin20 

p  po  o  2  o  p  y 


where  dots  indicate  derivations  with  respect  to  time. 

The  amplitudes  of  the  forced  osclillations*  calcula¬ 
ted  with  the  approximation  sin2©  »  20  can  be  read  in  Table  IV, 
for  the  same  satellites  considered  in  Table  III. 


-  33  - 


TABLE  IV  -  Amplitude  of  the  oscillations  of  the  satellite 
of  Pig. 1-5  (see  Eg. (1.7-2)  and  Table  III  . 


c Zi 


a 

(m) 

b 

(m) 

0 

(m) 

weight 

(hg) 

e  = 

amplitude 
,0-2  e  = 

10"1 

1 

1 

2 

500 

1° 

25' 

14° 

20' 

1 

1 

2 

1000 

1° 

25' 

14° 

20' 

1 

1 

2 

2000 

1° 

25' 

14° 

20' 

1 

1 

3 

1000 

45  s 

8° 

35' 

1 

1 

3 

2000 

45' 

8° 

35' 

2 

2 

4 

5000 

1° 

m 

CM 

14° 

20' 

2 

2 

4 

6500 

1° 

25' 

14° 

20' 

2 

2 

4 

8000 

1° 

25' 

14° 

20' 

2 

2 

4 

10000 

1° 

25' 

14° 

20' 

Table  IV  shows  that,  when  the  eccentricity  of  the 
orbit  is  not  very  small,  it  is  necessary  to  provide  the  sa¬ 
tellite  with  means  for  counteracting  the  natural  oscillations. 

R.P.  Haviland  (^)  states  that  a  synchronous  broad¬ 
cast  satellite  ,  weighing  500  lb  ,  with  a  power  demand  of 
1500  watts,  ''is  within  the  reach  of  present  technique' '  and 
that  11  a  single  50  kW  transmitter  (weight  =  10000  lb)  can  be 
installed  in  a  satellite  in  about  7  years ' '  . 

By  comparison  of  Table  I  and  Table  II  it  is  seen 
that  the  electromagnetic  torques  are  sufficient  to  counter¬ 
balance  the  natural  torques  on  a  satellite  of  this  kind. 
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The  damping  action  would  he  accomplished  practically, 
by  placing  a  turnstile  antenna  on  the  synchronous  satellite, 
with  the  axis  of  symmetry  of  the  antenna  normal  to  the  plane 
of  the  orbit  (i.e.  horizontal) , and  by  making  it  radiate  pov/er 
at  a  ratio  depending  sinusoidally  on  time.  In  order  to  obtain 
with  a  reasonable  expense  of  power  ,  torques  sufficient'  to 
overcome  the  disturbing  torques  indicated  by  (1.6-2)  and  ta  - 
ble  III  ,  a  long-arms  antenna  should  be  used  .  This  justifies, 
in  this  case  at  least,  the  assumption,  made  in  the  compilation 
of  table  I  ,  that  the  satellite  body  does  not  much  disturb  the 
field  radiated  by  the  antenna. 

The  device  described  above  would  allow  control 
about  a  single  axis  .  The  equatorial  position  of  the  satelli¬ 
te  would  suggest  the  possibility  of  obtaining  control  about 
the  other  two  axes  ,  by  simply  placing  a  big  permanent 
magnet  (  )  on  the  satellite  .  The  magnet t oriented  North-South 

(magnetic),  might  overcome  the  small  disturbing  torques  about 
the  other  axes. 

So  far  the  inertial  reaction  torques  that  a  satel¬ 
lite  must  overcome  ,  when  it  has  to*  keep  its  attitude  fixed 
with  i-espect  to  a  geocentric  reference  system,  have  been  eva¬ 
luated  for  the  particular  case  of  the  synchronous  satellite. 

Table  V  gives  the  torques  about  the  pitch  axis  , 

acting  on  satellites  orbiting  at  other  altitudes  .  The  moment 

2 

of  inertia  K  about  the  pitch  axis  is  assumed  to  be  100  kg  m  . 

P 

The  table  can  be  easily  extended  to  other  values  of  K  since 

3? 

the  torque  is  proportional  to  K 
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TABLE  V  -  Value  of  T  =  2Kew^  ,  for  different  altitudes  h 

o  o 

and  eccentricities  e  .  Tq  is  given  in  Joules  ,  h  in  km. 


h 

(km) 

CD 

ll 

0J 

ro 

To  «J) 

e  =  10”1 

e  =  2  10”1 

500 

2.4  10"6 

2.4  10“5 

ir\ 

1 

O 

CO 

1000 

2  10“ 6 

2  10” 5 

4  10”5 

2000 

1.4  10~6 

1.4  10” 5 

2.8  10"5 

5000 

VJ1 

o 

I 

5.4  10“6 

10.8  10”6 

10000 

1.8  10~7 

1.8  10"6 

-6 

3.6  10 

When  the  orbit  has  a  certain  inclination  with  re¬ 
spect  to  the  equator  ,  torques  about  the  roll  and  the  yaw 
axes  can  also  arise  in  a  geocentric  reference  system  ,  due  to 
the  regression  of  the  orbit  .  A  calculation  (**)  shows  that, 
even  for  satellites  of  large  moments  of  inertia,  these  torques 
have  values  of  a  few  ergs,  at  an  altitude  h  =  1000  km,  and  fall 
off  as  (R  +  h) ^  ,  (R  is  the  earth's  radius). 
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1 .8  -  Other  — — i  forces  acting  on  the  attitude  of  a, 

space  vehicle 

a)  Earth’s  magnetic  field. 

The  earth's  magnetic  field  can  influence  the 
attitude  of  a  space  vehicle  by  different  actions.  It  can 
interact  with  currents  induced  in  the  moving  body  ,  with 
currents  in  not  well  shielded  circuits  ,  with  induced 
magnetic  moments  and  with  permanent  magnetic  parts  of  the 
vehicle . 

As  it  generally  occurs  for  all  of  the  external 
actions,  the  values  of  the  torques  depend  largely  on  the 
particular  situation  considered  . 

One  can  have  an  idea  of  the  order  of  magnitude  of 
the  magnetic  actions  by  considering  some  practical  cases. 

The  attitude  of  the  satellite  Tiros  I  was  stabi¬ 
lized  by  means  of  a  rotation  about  the  axis  of  maximum 
moment  of  inertia  .  It  was  noticed  that  the  direction  of 

i  9 

this  axis  changed  with  time.  This  motion  was  explained  (  * ) 
by  taking  into  account  the  interaction  of  the  magnetic 
parts  of  the  satellite  with  the  earth's  magnetic  field  . 

The  calculated  magnetic  moment,  directed  along  the  spin 
axis  ,  was  896  ergs/G  .  At  the  altitudes  of  Tiros  I  the 

magnetic  field  is  about  0.4  G  y  thus  the  resulting  va¬ 
lues  of  the  torque  turn  out  to  be  of  the  order  of  hundreds 
of  ergs  . 

For  the  Explorer  IV  ,  the  calculations  of  G.  Co¬ 
lombo  ,  based  on  the  average  diurnal  variation  of  the  an- 
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gular  momentum,  of  the  satellite  ,  as  sume^a r  jiaximum  magnetic 
torque  of  300  ergs  ;  .  He  also  reports  that  the  amount 

of  magnetic  torque  predicted  for  the  Orbital  Astronomical 
Observatory  ,  is  many  hundred  ergs. 

These  values  of  the  magnetic  torques  seem  high  , 
when  compared  with  the  electromagnetic  torques  of  Table  I  . 

It  must  be  observed  however  ,  that  the  magnitude  of  the  earth's 
magnetic  field  falls  off  with  the  third  power  of  the  distance 
from  the  center  of  the  earth  .  Thus,  satellites  at  altitudes 
higher  than  those  of  the  Explorer  IV  (minimum  altitude  less 
than  300  km.), of  the  Tiros  I  (altitude~700  km  )  and  of  the 
O.A.O.  (altitude<“800  km)  ,  can  be  subjected  to  magnetic  tor¬ 
ques  weaker  than  those  reported  above  .  Besides,  the  inter¬ 
planetary  magnetic  field  is  very  weak  (')  ,  and  the  magnetic 
torques  acting  on  vehicles  travelling  far  from  the  earth  and 
other  planets  would  also  be  very  small. 

An  indirect  action  of  the  magnetic  field  on  the 
attitude  of  a  satellite  is  also  to  be  considered  :  a  conducting 
body  rotating  in  a  magnetic  field  is  slowed  down  by  the  interac¬ 
tion  of  the  magnetic  field  with  the  induced  currents. 

For  this  reason,  spin  stabilized  satellites  lose 


( ' )  Due  to  emission  of  plasma  clouds  or  of  coronal  gas  from 

j 

the  sun,  the  solar  magnetic  field  is  transported  into  the 

interplanetary  space, giving  rise  to  .“magnetic  fields  of 

—  3  -5 

10  -10  G  in  regions  near  the  sun.  Beyond  the  reach  of  the 

solar  corpuscolar  stream  a  galactic  field  of  the  order  of 
10“ ?G  is  expected  to  prevail  through  outer  space  (  ) 
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rotation  speed  and  consequently  stability;  of  the  .spin  axis. 

A  study  of  this  braking  action  on  vehicles  of.  simple 

22 

geometrical  shape  has  been  made  by  Wilson  .(  )  . 

The  slowing  torque  depends  on 'the  electrical  proper¬ 
ties  of  the  material  constituting  the  satellite  and  on. its  shape. 
Pox'1  spherical  or  cylindrical  satellites  it  is  proportional  to 
the  spin  velocity  .  Such  an  effect  was  noted  for  the  Vanguards  I 
and  II  and  for  the  Tiros  I  . 

The  spin  velocity  u>  of  the  Vanguard  I  (satellite  of 

s 

spherical  shape  ,  diameter. =  16  cm),  was  found  to  vary  with 
time,  according  to  the  law  (')  : 

u)g  =  2.72  exp[-5.03  10"8(t-tQ)]  r/e 

where  t  is  expressed  in  seconds,  and  t  denotes  the  time  at 
which  the  satellite  was  placed  in  orbit. 

The  moment  of  inertia  I  of  the  satellite  was 
6.9  ICT3  kg  m2. 

Hence  the  slowing  torque  C  results  to  be 

C  -  3.48  10"10u>  J/s 

c 

(The  proportionality  to  wg .indicates  the  .action  of  the  earth's 
magnetic  field) .  •  .  .  • 


(')  The  data  for  the  Vanguard  satellites  is  taken  from  Wilson's 
paper  (see  Ref. 22) 
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Thus  ,  at  the  beginning  of  its  life,  the  satellite 

was  submitted  to  a  slowing  torque  of  9.5  10-1®  J. 

A  more  important  case  is  the  Vanguard  II  ,  a  spheri¬ 
cal  satellite  of  9.75  kg  with  a  diameter  of  0.51  m  .  A  well 

determined  angular  speed  was  important  not  only  for  the  sta¬ 
bilization  but  also  for  the  better  performance  of  the  photogra¬ 
phic  apparatus. 

Its  spin  velocity  was  found  to  vary  with  the  law 


The  moment  of  inertia  about  the  spin  axis  was 
1 .977  kg  m^  .  Hence  : 

C  =  0.3205  10“7  u  J/s 

s 

Then  the  initial  slowing  torque  acting  on  the  satellite  was: 

0*8  10“9  J 

o 

When  an  external  torque  is  applied  to  a  spinning 

satellite  ,  it  gives  rise  to  a  precession  of  the  spin  axis  . 

The  magnitude  of  the  precession  is  proportional  to  the  torque, 

to  the  duration  of  (the  action  >  and  inversely  proportional  to 

23 

the  angular  momentum  of  the  satellite  about  the  spin  axis  (  ) . 

This  latter  quantity  may  be  taken  as  an’  index  of  the  stability 
of  the  satellite  and  in-  the  case  of  the  two  Vanguards  was 
reduced  to  one  half  of  its  initial  value  within  160  and  50  days 
respectively  . 
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Wilson' 3  analysis  shows  that  the  slowing  torque  and 
the  moment  of  inertia  ,  for  spherical  bodies,  are  proportional 
to  the  same  power  (the  fifth)  of  the  vehicle's  radius  .  Thus 
for  satellites  of  larger  dimensions  hut  with  similar  distri¬ 
bution  of  mass  and  materials  ,  as  the  Vanguards,  the  time  in 
which  the  spin  velocity,,  or  the  inherent  stability  ,  is  reduced 
to  one  half  ,  is  of  the  same  order  of  magnitude  (at  the  same 
altitudes)  . 

The  slowing  down  of  the  spin  velocity  was  taken 
into  account  in  the  post-launch  performance  of  the  Tiros 
satellites . 

The  Tiros  I  (of  cylindrical  shape  ,  diameter  =  1.05 
weight  5  kg)  began  its  working  life  with  a  spin  velocity  of 
12  r/min  .  After  5b  days  this  value  was  reduced  to  9.4  r/min. 

A  pair  of  rockets  ,  placed  around  the  satellite,  brought  up 
the  rotation  velocity  to  12.85  r/min.  Two  other  pairs  of 
rockets  were  necessary  during  the  life  of  the  satellite. 

If  one  assumes  that  also  in  this  case  u  varies 

8 

exponentially  with  time 


U) 

s 


2 

and  takes  into  account  the  value  12.4  kg  m  of  the  moment  of 
inertia  about  the  spin  axis  ,  one  finds 


G  =  8  10“7  J 

o 


for  the  value  of  the  initial  torque 


b)  Solar  radiation  pressure 


The  3olar  radiation  pressure  depends  on  the  angle  of 
incidence  of  the  radiation  and  on  the  reflectivity  of  the  sur¬ 
face  of  the  vehicle.  For  complete  reflection  the  pressure  can 
reach  a  value  of  9  10^  dyne/cm^  at  the  earth's  distance  from  the 
sun  .  For  complete  absorption  ,  which  is  generally  considered 

a  good  approximation  ,  at  this  distance  ,  the  sun  pressure  P 

-4/2  3 

reaches  a  value  Pge4.5  10  dyne/cm  . 

In  the  case  of  complete  absorption  the  torque  Tg 

which  can  result  is  the  product  of  : 

the  area  A  presented  to  the  sun, 
the  radiation  pressure  Pg 

and  the  distance  D  between  the  center  of  pressure  of 

the  area  A  and  the  projection  of  the  center  of  mass 

of  the  vehicle  on  A. 

With  the  above  value  of  P  , 

s 


T  =4.5.10-^  AD  Joules 

s 

2 

where  A  is  in  m  ,  and  D  in  m. 

It  is  seen  that  for  a  nearly  symmetrical  vehicle  the 
perturbing  action  of  the  solar  pressure  on  the  attitude  can  have 
very  small  values . 
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c)  Meteorites  impact. 

The  impact  of  a  meteorite  on  a  space  vehicle  can 
cause  an  impulsive  change  of  the  angular  momentum  of  the 
vehicle  .  Whipple  (2^)  gives  the  estimated  distribution  of 
the  masses  and  the  velocities  of  small  meteorites  near  the 
earth  .  He  also  gives  the  calculated  mean  number  per  day 
of  meteorites  striking  a  3  m  diameter  sphere  ,  placed  in  the 
space  near  the  earth  .  This  number  is  derived  with  the  assump¬ 
tion  of  50  %  shielding  by  the  earth. 

Let  us  consider  a  spherical  space  vehicle  of  this  size. 
The  second  column  of  table  VI  gives  the  number  of  particles, of 
masses  in  the  range  indicated  in  the  first  column,  which  strike 
the  vehicle  per  day  .  The  values  are  derived  from  Whipple's 
data  ,  by  neglecting  the  shielding  effect  of  the  earth. 

The  third  column  gives  the  angular  momentum  I, 
absorbed  by  the  vehicle  when  a  complete  unelastic  impact  and 
a  distance  d  =  1  m  between  the  trajectory  of  the  meteorite 
and  the  center  of  mass  of  the  vehicle  are  assumed  .  The  rela¬ 
tive  velocity  of  the  meteorite  with  respect  to  the  vehicle 
has  been  assumed  to  have  the  value  of  80  km/s  (near  to  the 
estimated  maximum)  . 

For  the  mass  of  the  meteorite  we  used  the  average 
between  the  figures  given  in  each  row. 
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TABLE  VI  -  Data  on  meteorites  impact  (see  text) 


range  of  masses 

(gr) 

number  per 
day 

I=change  of  an¬ 
gular  momentum 
(j-sec/unit  of  d) 

0.25  -  9.95  10" 2 

4.3  10“5 

14 

9.95  l0_2-3.96  10“2 

9.8  10“ 5 

4.8 

3.96  10"2-1.58  10“2 

2.5  10“4 

2.2 

1.58  10"2-6.28  10“3 

6.2  10“4 

9.1  10"1 

6.28  10“  3-2. 5  10“  3 

1.6  1 0” 3 

3.5  io”1 

2.5  10"3-9 .95  10"4 

3-9  10“ 3 

1,  .4  10-1 

9.95  10"4-3.96  10“4 

9.8  10“3 

4.8  10“2 

3,96  10-4-l .58  10"4 

2.5  10“2 

2.2  10“ 2 

1.58  10"4-6.28  10“5 

6.2  10"2 

9.1  10"3 

6.28  10“5-2.5  10“ 5 

1.6  10"2 

3.5  10“3 

The  table  can  be  extended  to  other  sizes  of  the 
vehicle,  by  considering  that  the  number  of  particles  striking 
a  sphere  is  proportional  to  the  square  of  the  diameter. 
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It  is  seen  from,  table  VI  that,  especially  for  satel¬ 
lites  of  long  life  (')  ,the  probability  of  impact  with  meteo¬ 
rites  of  considerable  mass  cannot  be  excluded  .  The  magnitude 
of  the  change  of  angular  momentum  can  be  large  .  Thus  an 
attitude  control  device  should  not  ignore  the  effect  of 
meteorites  impact  ,  and  should  be  able  to  exert  .  large  torques, 
especially  if  this  perturbing  action  has  to  be  balanced  in  a 
short  time . 

Por  example,  let  us  consider  a  spherical  vehicle, 
with  diameter  —  1  m,  which  ,  initially,  has  its  attitude  fixed 
with  respect  to  an  inertial  reference  frame. 

-4 

Let  us  assume  that  a  meteorite  ,  having  mass  =  7.10 
grams  and  travelling  with  a  velocity  =  80  km/s  relative 
to  the  vehicle  ,  impinges  onto  the  vehicle  .  Prom  table  VI 
it  is  seen  that  the  average  number  of  impacts  with  a  similar 
meteorite  is  one  every  2-3  years,  for  a  vehicle  of  diameter 
=  1  m.  Let  us  also  assume  that  the  distance  d  between  the 
center  of  mass  of  the  vehicle  and  the  trajectory  of  the  me¬ 
teorite  is  0.25  m  ,  and  that  the  impact  is  unelastic  . 

Then  the  vehicle  ,  receiving  an  impulsive  acceleration,  begins 
to  rotate  about  an  axis  normal  to  the  direction  of  the  meteori¬ 
te  ,  with  an  angular  velocity  m  =  g  ,  where  K  is  the  moment 
of  inertia  about  this  axis  .  In  this  case  ,  I  *  1.2  10  a  s 
(see  table  VI  column  III)  . 

If  a  device  must  stop  this  motion  before  the  attitu¬ 
de  deviation  reaches  a  certain  angle  a,  the  control  torque  T, 
assumed  to  be  constant  during  the  control  action  ,  must  be  : 


(')  Por  instance,  it  is  assumed  that  an  active  communications 
satellite  would  be  economically  convenient,  only  if  its  life 
could  be  some  yers  long. 
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TC  2  fa 

For  the  example  discussed  ,  table  VII  gives  the 

values  of  T  ,  for  several  values  of  a  and  K  .  The  table 
c 

also  shows  the  time  t  ,  necessary  to  stop  the  rotation  of 

s 

the  vehicle,  with  the  indicated  control  torques  . 

TABLE  VII  -  Values  of  T 

c 


a 

II  =  100  kg  m2 

T  (J)  t  (sec) 

0  8 

2 

K  =  1 000  kg  m 

T  (J)  t  (sec) 

0  S 

5° 

0.8  10"5 

1.5  1 0'3 

0.8  10“6 

1.5  104 

1° 

4.13  10“ 5 

2.9  102 

4.13  10"'6 

2.9  103 

10' 

2.48  10~4 

1 .93  103 

2.48  10"5 

1.93  104 

1  1 

2.48  10“3 

1.93  102 

-4 

2.48  10 

1.93  103 

The  table  is  easily  extended  to  other  values  of 

K  and  T  . 

After  the  motion  of  the  vehicle  has  been  stopped, 
the  control  device  has  to  readjust  the  attitude  of  the 
vehicle  ,  by  making'  it  rotate  of  an  angle  a  in  the  opposite 
direction.  To  have  an  idea  of  the  time  t  necessary  to  ob¬ 
taining  the  rotation  with  a  given  controlling  torque  T  ,  let 
us  consider  the  following  procedure.  The  device  exerts  the 
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torque  Tc  for  a  time  tQ.  Then,  for  an  equal  time,  the  device  keeps 


acting  with  the  same  but  opposite  torque.  At  the  end  of  these  two 
phases  the  vehicle  has  accomplished  a  rotation  a  and  has  been 
stopped.  The  angle  a  is  given  by: 


1  Tc  .2 

j  r  *t 


where  t  =  2t  In  Pig. 1.6  the  angle  a  (measured:  in  jcadiants)  is  plotted 
r  o 


U  (nad ) 

(Tc  -•fboocrg)  (T,;=fOOcr3)(Tc=^erg) 


versus  t  (  in  seconds)  for  different  values  of  K  , 

The  three  scales  of  the  ordinate  axis  correspond  to 

T  =  106  ,  T  =  10" 5  ,  T  =  10~4  j 
c  c  ’  c  * 

Of  course  ,  during  the  oime  t  ,  the  device  must 

r 

exert  a  supplementary  torque,  in  order  to  balance  the  actions 
of  other  sources  of  disturbance. 

d)  Atmospheric  drag. 

The  motion  of  a  satellite  in  the  earth  atmosphere 

gives  rise  to  an  aerodynamic  pressure  acting  on  the  vehicle. 

If  the  centers  of  pressure  and  of  mass  of  the  vehicle  are  not 

in  the  same  line  with  the  direction  of  the  orbital  motion  ,  a 

torque  results  applied  to  the  vehicle  •  The  aerodynamic  pressure 

P,  on  an  orbiting  vehicle  varies  very  largely  with  the  altifcu- 
u  25 

de  .  Prye  (  )  indicates  the  following  values  of  P,  : 

d 

at  an  altitude  of  1 30  km  P  =  4.5  104  dyne/m^ 

®  _  'J  p 

at  an  altitude  of  560  km  .  P  =  1  .4  10  dyne/m 

d 

It  is  seen  that  the  aerodynamic  pressure  can  play 
an  important  role  with  regard  to  the  attitude  of  a  vehicle, 
at  low  altitudes.  At  high  altitudes  however  the  aerodynamic 
torques  can  be  neglected  in  comparison  with  other  causes 
affecting  the  attitude  .  For  example  it  is  to  be  recalled  that 
the  solar  radiation  pressure  is  of  the  order  of  4  -  8  dyne/m  . 

e)  Motion  of  internal  parts. 

The  motion  of  internal  parts  can  influence  the  at- 
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titude  of  a  space  vehicle  * 


For  example  the  change  of  the  position  of  a  part  of  a 
vehicle  with  respect  to  the  main  body  ,  or  the  acceleration 
of  rotating  parts,  can  induce  a  torque  on  the  rest  of  the 
vehicle  through  inertial  reaction. 

Gyroscopic  torques  can  arise  when  the  attitude  of  a 
vehicle,  containing  rotating  parts  ,  has  to  he  changed  (or  must 
stay  fixed  in  a  non  inertial  reference  frame)  =, 

To  have  some  numerical  examples  ,  let  us  consider 
the  cases  :  1)  a  vehicle,  whose  only  moving  parts  are  the 
electrical  circuits  switches  ,  and  2)  a  vehicle  containing 
an  electric  turbo-engine  . 

In  the  first  case  one  can  see  that  the  disturbing 
actions  induced  hy  the  the  moving  parts  are  very  small.  In 
fact,  let  us  consider  a  small  mass  m  moving  along  an  arc  s  , 
of  radius  r,  about  the  center  of  mass  of  the  mainbody  . 

The  vehicle  ,  by  reaction  ,  will  make  a  rotation 
a  in  the  opposite  direction  .  If  K  is  the  moment  of  inertia 
of  the  vehicle  about  the  axis  of  rotation  ,  a  is  given  bys 

smr 

a  =  -p. 

For  instance  ,  if  :  s  =  1  cm  ,  m  =  1  g  , 

2 

r  =  1  m  ,  K  =  TO  kg  m  ,  a  turns  out  to  be  about  0»2  se¬ 
conds  of  arc. 

For  the  second  case,  let  us  consider  a  vehicle 

carrying  a  turbo-electric  generator  ,  with  a  rotor  having  moment 

2 

of  inertia  =1  kg  m  ,  rotating  at  20.000  r/min  . 


Let  us  suppose  that  the  attitude  of  the  vehicle  must 

he  changed  with  respect  to  an  inertial  reference  frame:  the 

vehicle  has  to  rotate  about  an  axis  making  an  angle  a  with 

the  axis  of  the  rotor.  The  angular  velocity  of  this  rotation 

is  oj  .  In  this  case  a  gyroscopic  torque  T  results  applied  to  the 
v  g 

vehicle. 


T  =  u  H  sina 
S  v 

where  H  is  the  magnitude  of  the  singular  momentum  of  the  rotor 

with  respect  to  its  axis. 

—  3  —1 

If  =  10  sec  (case  of  a  body  orbiting  around 
the  earth  at  an  altitude  h  =  1000  km,  with  an  axis  constantly 
directed  towards  the  earth) ,  and  a  «  1°  ,  T  turns  out  to  be 

_3  £ 

equal  to  ~  6 . 1 0  J  . 

The  two  preceding  examples  shows  that  the  motion  of 
internal  parts  can  give  rise  to  torques  applied  to  a  vehicle, 
which  vary  in  a  very  wide  range,  according  to  the  particular 
situation. 

1.9-  Conclusions  on  the  applicability  of  the  electromagnetic 
torques. 

In  sections  1.6  ,  1.7  said  1.8  the  order  of  ma¬ 
gnitude  of  the  torques  applied  to  a  space  vehicle,  due  to 
external  forces  or  to  .the  motion  of  internal  parts,  were 
discussed. 
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In  general  the  magnitudes  of  the  disturbing  actions 
depend  largely  on  the  particular  situation,  namely  on  the  dimen¬ 
sions  of  the  vehicle,  on  its  distribution  of  masses,  on  its 
shape,  and  on  the  characteristic  of  the  missions  to  be  accom¬ 
plished.  In  particular,  for  an  earth's  satellite  the  magnitudes 
of  the  forces  effecting  its  attitude  vary  much  with  the  altitude 
of  the  orbit. 

Especially  when  the  altitude  is  high  enough,  so  that 
the  influence  of  the  gravitational  field, of  the  magnetic  field 
and  of  the  atmosphere  is  not  of  great  importance  ^see  sec. 1.6  a) 
and  1.8  a),b)^,  the  torques  which  can  be  obtained  by  means  of 
emission  of  circularly  polarized  radiation  are  sufficient  for 
balancing  external  forces.  It  seems  therefore,  that  the  elec¬ 
tromagnetic  torques  can  be  employed  for  the  control  of  the 
attitude  of  a  stabilized  vehicle.  One  should  take  int:o  parti¬ 
cular  account  the  possibility  of  impact  with  the  meteorites.  As 
shown  in  section  1.8  b),  the  impact  wijth  a  meteorite  of  sub¬ 
stantial  mass  can  be  considered  a  very  rare  event.  However,  in 
some  cases  of  long  life  vehicles,  it  may  be  necessary  to  fur¬ 
nish  the  vehicle  with  a  control  device  to  counteract  this 
occurrence.  Sec. 1.7  b)  3ho\vs  that  relatively  high  changes.of 
angular  momentum  of  the  vehicle  can  ensue. -The  torques  to  be 
exerted  by  the  controlling  device  depend  on  the  tolerated  at¬ 
titude  deviation  a,  and  on  the  permitted  length  of  the  time  t 
wjthin  which  the  control  device  must  readjust  the  attitude. 

The  possibility  of  employing  electromagnetic  torques,  in  this  case 
depends  on  the  values  of  a  and  tj.,  as  well  as  on  the  power  and 
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wavelength  available. 

In  section  1 .7  the  particular  case  of  a  synchronous 
satellite  was  considered,  together  with  a  possible  scheme  of 
a  single  axis  control.  For  a  three  axis  control  three  turnstile 
antennas  should  be  placed  on  the  satellite,  each  one  control¬ 
ling  the  attitude  about  a  certain  axis  . 

More  simply,  a  three  axis  control  could  be  achieved 
by  means  of  a  set  of  three  half  wave  antennas,  disposed  as  in 
Figure  1-7. Bach  of  the  three  possible  pairs  of  antennas,  can 


Fig.1-<7  -  Antennas  for  a  three  axis  control. 
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form  a  turnstile  antenna.  Different  frequencies  could  be  used 
for  the  three  turnstile  antennas,  in  order  to  permit  an  easy 
separate  feeding. 

Other  types  of  antennas  could  be  considered.  However 
it  is  to  be  expected  that  in  general  the  reaction  torques  ac¬ 
ting  on  antennas  radiating  with  circular  polarization  would  be 
proportional  to  the  radiated  power  and  to  the  wavelength.  Thus 
it  seems  that  helical  antennas  or  spiral  antennas,  whose  di¬ 
mensions  should  be  of  the  order  of  the  wavelength,  in  order  to 
have'  a  good  electrical  efficiency,  would  be  more  massive  and 
cumbersome  than  the  simple  turnstile  antenna. 

Fig.W6  shows  in  what  limits  a  device,  capable  of 
supplying  weak  torques,  can  be  employed  when  the  attitude  of  a 
space  vehicle  must  be  changed.  In  fact  it  shows,  for  different 
values  of  the  torques  available  and  of  the  moment  of  inertia 
of  the  vehicle  about  the  rotation  axis,  the  time  necessary  to 
obtain  a  given  rotation  a.  Of  coursjs,  during  the  period  of  the 
attitude  correction,  in  addition  to  the  torque  necessary  to 
change  the  attitude  of  the  vehicle,  considered  free  from  ex¬ 
ternal  actions,  the  device  has  to  supply  an  amount  of  torque 
sufficient  to  counteract  the  external  actions,  arising  from 
the  various  causes  discussed  in  sections  1.6  and  1.8. 


1.10  -  The  electromagnetic  torque  as  a  disturbing  action. 

The  different  sources  of  disturbing  torques,  listed 
in  sections  1.6  and  1.8,  as  well  as  other  of  minor  importance. 
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have  been  widely  studied  in  the  literature.  So  far,  the  torques 
which  can  arise  due  to  the  emission  of  circularly  polarized 
radiation  seem  not  to  have  been  considered. 

The  power  transmitted  by  the  antennas  placed  on  the 
satellites,  launched  up  to  now,  was  sufficiently  low,  so  that 
the  torques  arising  in  the  case  of  circular  polarization  were 
very  weak. 


As  an  example,  let  us  consider  the  case  of  the  satel¬ 
lite  Explorer  VII,  which  carried  a  turnstile  antenna  radiating 
600  mW  at  20  Me. 

.  .  .9 

By  applying  Eq.(1.4-13)  a  torque  H  5  10  J  turns 

out  to  be  applied  to  the  satellite. 

Explorer  VII  was  stabilized  by  means  of  a  rotation 
about  the  axis  normal  to  the  plane  of  the  turnstile  antenna. 

For  a  comparison  with  the  electromagnetic  torque  one  could  con¬ 
sider  the  torque  T  which  was  actually  applied  to  the  satellite 
about  this  axis.  This  can  be  derived  from  the  slowing  down  of 
the  spin  velocity. 

By  assuming  an  exponential  decrease  of  the  spin  velo¬ 


city  u  : 
s 

one  has: 


-at 

u>  =  to  e 
s  0 


Where  K  is  the  moment  of  inertia.  An  approximate  value  of  K  is 
obtained  by  assuming  the  satellite  to  be  of  spherical  shape, 
with  a  homogeneous  distribution  of  mass.  Then,  using  the  data 
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derived  from.  (^°),  one  has: 

2 

K  =  ~  2.4  kg  m 

the  initial  spin  velocity  was  3»0  r/m.  After  50  days  it  was 
reduced  to  350  r/m.  One  obtains  for  Mq  (initial  torque): 

M  =  ~  6  10  '  J. 

0 

This  is  two  orders  of  magnitude  greater  that  the  electromagnetic 
torque . 

In  the  future  however  high  levels  of  power  will  be 
necessary,  especially  for  communications  satellites  (as  indi¬ 
cated  in  section  5). 

27 

It  has  been  pointed  out  (  )  that  for  earth-space 

communications:  * 'if  we  have  to  avoid  serious  Faraday  fading 
for  a  wide  range  of  signal  frequencies  circular  polarization 
or  special  reception  technique  are  both  desiderable  and  neces¬ 
sary  . ' ' . 

In  sections  1.6  and  1.8  it  was  shown  that  the  torques 
applied  by  reaction  to  a  space  vehicle,  when  circularly  polari¬ 
zed  waves  of  high  power  are  radiated,  can  be  of  the  same  order 
of  magnitude  as  the  other  disturbing  torques  applied  to  the 
vehicle.  Therefore  in  the  prediction  of  the  post^launoh  perfor¬ 
mance  of  a  satellite  with  an  antenna  radiating  circular  pola¬ 
rization,  the  action  of  the  electromagnetic  torques  on  the 
attitude  of  the  vehicle  will  have  to  be  taken  into  account. 
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2  -  TORQUE  ON  A  RECEIVING  ANTENNA 


2.1  -  introduction 


In  the  following  sections  a  possible  second  way  will 
be  discussed  of  applying  the  properties  of  elliptically  pola¬ 
rised  electromagnetic  waves,  in  order  to  exert  a  torque  on  a 
space  vehicle  ,  to  control  its  attitude. 

let  us  consider  the  simple  case  of  elliptically  po¬ 
larized  radiation  impinging  on  an  antenna,  placed  in  the  free 
space  and  far  from  other  material  bodies  .  The  angular  momentum 
transferred  from  the  waves  to  the  antenna  per  unit  time  (that 
is  the  torque  exerted  by  the  radiation)  will  be  evaluated. 

The  torque  is  equal  and  opposite  to  the  outgoing  flux  of  the 
angular  momentum  of  the  radiation  through  the  surface  I  of 
a  sphere  (')  ,  centered  at  the  antenna  .  The  flux  will  be 
evaluated  by  taking  into  account  both  the  incident  and  re-ra¬ 
diated  waves. 

To  get  an  idea  of  the  order  of  magnitude  of  the  tor¬ 
que  which  could  be  obtained  in  this  way,  let  as  consider  a 
turnstile  antenna  under  the  action  of  a  plane  circularly  pola¬ 
rised  wave,  propagating  in  a  direction  perpendicular  to  the 
plane  of  the  antenna  .  The  antenna  is  assumed  to  be  resonant 


(') 


see  sec.  1.2 
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at  the  wavelength  X  of  the  incident  wave  .  Each  component 
antenna  has  a  length  =  and  its  terminals  are  closed  through 
a  resistance  E  equal  to  the  radiation  resistance. 

With  reference  to  Pig. 2-1,  let  the  plane  (x,  y)  he 
the  plane  of  the  turnstile  antenna,  0  the  origin  of  the  axes, 
coincident  with  the  center  of  the  antenna.  The  two  half-wave 
component  antennas  lie  on  the  x  and  y  axes  respectively. 


Pig. 2.1 
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Here  and  in  the  sequel  all  elements  of  the  y  antenna 
will  be  labelled  by  the  subscript  1,  and  all  elements  of  the  x 
antenna  by  the  subscript  2. 

let  the  electric  and  magnetic  fields  of  the  incident 
circularly  polarized  plane  wave  be  given  by: 


(2.1-1) 


E 

-o 


E0  (i  +  W 


i(wt  -  kz) 
e 


(2.1-2)  Hj  -  J  kA^  -  ^2  U  -  i  i)  ei(ut  "  kz) 

Here  i  denotes  the  imaginary  unit,  Z  the  characteristic  impedan¬ 
ce  of  the  free  space,  and  k  =  2ti/X  the  wave  number. 

The  mean  value  of  the  component  on  the  z-axis  of  the 
torque  acting  on  the  antenna  is  given  by  the  expression  ( •■) : 


(2.1-3) 


M  =  *  Re 


(EAn.k) (D*.n)  +  (H^n.k) (B*.n) 


r  d£ 


where  I'  represents  a  spherical  surface  with  radius  r  centered 
at  0.  E,  H,  D  and  B  are  the  field  vectors;  n  is  the  outward 
normal  to  I  . 


(')  see  sec. 1.2  Eq.( 1.2-7) 
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2.2  -  Torque  on  a  turnstile  ant enna  under  the  action  of  a  cir¬ 
cularly  -polarized  plane  wave . 

The  incident  wave  gives  rise  to  currents  on  the  two 
half-wave  component  antennas. If  the  branches  of  the  antenna  are 
assumed  to  be  infinitely  thin,  the  expressions  of  the  currents 
are: 

(2.2-1)  Ig  =  I  e^^cos  (kx ) 

(2.2-2)  1^  =  Ho  eiwtcos(ky) 

The  total  electromagnetic  field  E,  H  is  the  sum  of 
the  incident  field  plus  the  field  re-radiated  by  the  two  dipo¬ 
les.  Thus 


(2.2-3) 


(2.2-4)  H  =  H  +  H  i  +  H_  i 

-  “0  Icp-xpl  2q>  (p2 


where  n,  irt.  ,  i  ,  and  n,  i_„,  i  0  denote  the  unit  vectors  of  two 

— ’  =01  ’  — cp  1  —  -Q2  — cp 2 

spherical  coordinate  systems,  having  the  polar  axes  on  the  x  and 
y  axis  respectively  (see  Pig. 1-2), 
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(2.2-5) 


(2.2-6) 


(2.2-7) 


(2.2-8) 


(2.2-9) 


(2.2-10) 


(2.1-3), 


ZI  cos 
0 


\  oos«1 


1 0 


"Sir  ”  “ilnte "  e 


i(wt  -  kr) 


E 


1r 


ZI  X 
o 

x - Ti  3in 


|  cosOt 


i(wt  -  kr) 


H 

Icp 


I  C0S|5  oosfl  .  /  .  .  X 

0 _ \2 _ y_  x(wt  -  kr) 

?7ir  sln5“ 


ZI  cos  Us  cosO 

E  _  j_  o _ \2 _ 2/ 

20  27ir  sin5 


71 


i(ut  -  kr) 


_  ZIoX  .In  a\  i(u>t  -  kr) 
E2r  *  '”5  sln(?  00s®2 )  8 


71 


I  COS  •*  cos©,.  /  . ,  . 

ti  =  ,  _o _ \2 _ 27  i(wt 

2cp~  2iitr  ™  sinS^ 


-  kr) 


By  substituting  Eq_s.  (2.2-3)  and  (2.2-4)  into 
one  derives 
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(2.2-11)  H  = 


e  r 

Re 


'-•o+El0i<51+E20i02^A-*- 


E*  .n+E*  +E' 
-o - 1r  2r 


vr 


!(H  +H  i  ,+H  i  _ 
L  ~o  Icp-Tpl  2cp — cp 2 


)/\n.kj 


(H^.n)dl 


Let  us  now  evaluate  the  expression  (2.2-11).  From 
Fig. 2.1  the  follovdng  relations  can  be  derived: 

(2.2-12)  a  =  sinO^  coscp^  k  +  sin©^  sincp1  i  +  cosO^  ^ 

(2.2—13)  n  =  sin©2  coscp2  k  -  sin©^  sincp2  £  +  cos©2  i 
Hence  one  haa,  for  the  products  appearing  in  Eq. ( 2.2-1 1 ) , 


(2.2-H)  i/i  n.k  =  cos©1  =  -  sin©2  sirup  2 

(2.2-15)  i^n.k  =  -  sin©^  sincp1  *  -  coa©2 
(2.2-16)  i.n  =  sin©^  sincp^  =  cos©2 
(2.2-17)  j..n  =  cos©^  =  -  sin©2  sincp2 
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d£t 


(2.2-18)  (i^A  n)  .k  =  i^ .  k  =  cos«1  coscp., 


(2.2-19)  (i?p2A-^*-  =  i©2*-  =  C0S°2  C0S<P2 


(2.2-20)  i^^n.k  =  3incp1 


(2.2-21)  i^gAfi.k  =  sincp2 


By  using  E^.(2.1-l)  and  (2.1-2)  and  taking  into  account 
Eq^s .(2.2-14)  to  (2.2-17),  one  has: 


(2.2-22) 


Re 


(E^A^.k) 


( cosO^ 


sinO^ 


sirup  .j  -  cosO^  sinO^  sirup^)  »  0 


Analogously: 


Re 


/H  a  n.k)  B*.n 
\— o  —  —I  —  o  — 


=  0 


Eq.. ( 2.2-11 ),  by  the  help  of  Eqs.( 2.2-18)  through  (2.2-21),  beco¬ 
mes  .* 
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e  r  \  j 
o  -  ) 


(2.2-23)  Mz  =  Ee^J(B1^siap1  +  E^sincp^.n  dl  + 


+  (^.n.k)(E*r  +  E*r)dl  +  l(El0sincp1  +  E^siiupg)  (E*r  +  E^r)d£ 


Hr  j 

+  Re  {  (H,  cosO  cosm  +  H_  cos©„cos(p0)  h*.  n 

d  j  ?  » cp  1  1  <:cp  d  d  0  — 


Let  us  consider  the  integral,  appearing  in  Eq^(2.2-23) : 


(2.2-24)  A^=  J[EoAn.k](E^r+  E^dl'  =  r  dq^  /  [l^n.kjB^sin^dfi^ 


2u  *.71 


'0  v0 


*2n  ,Kit 

,f  .  i 


+  r2/dCp2|  [^oA-*-]^2r  3in°2  d0; 


0  vo 


From.  (2.1-1),  (2.2-14)  and  (2.2-15)  one  derives! 


(2.2-25)  E^n.k  =  Eq  ei^u,"b  “  'tr'sin®1C0st^1  ^cosO^  -  i  sinO.^  sirup  .^1 
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(2.2-26)  EQAn.k 


E  ei(ut  -  krsin©2cos<p2) 
o 


ainOgSincpg  -  i  cos©2 


The  second  term  of  the  sum  on  the  right  side  of  (2.2-25) 
is  a  symmetric  function  with  respect  to  © 1  =  n/2.  This  function 
multiplied  by  E^sinO^  gives  rise  to  an  antisymmetric  function, 
whose  integral  over  ©^  between  0  and  n,  in  Eq.(2.2~24),  vanishes. 
The  same  result  is  obtained  by  considering  the  first  term  of  the 
sum  on  the  right  side  of  (2.2-26). 

By  substituting  (2.2-25)  and  (2.2-26)  into  (2.2-24),  by 
taking  into  account  (2.2-6)  and  (2.2-9),  one  obtains: 


E  ZI  Xe 
o  o 


ikr 


(2.2-27)  A^=  -i  --g-- - |  sin^cos©]  sin©  cos©  d©  j  e 


2 n 

-ikr  sin©  cos(p 


The  integral  over  cp  is  equal  to  2n  JQ(kr  sin©).  Then 


(2.2-28) 


e  r 

-rHe  Ar 


rE  ZI  X 


.  i’Ji  4JX  A  /  \ 

- — g — --sin(kr)/  sin^cos©j 


JQ(kr  sin©) 


sin©  cos© 


With  the  substitution: 


(2.2-29) 


the  integral  in  (2.2-28) 


becomes: 


r 

(2.2-30)  |  sin 

\) 

GO 


injj;  cos«J  Jq 


(hr  ain©)sinfl  cosQ  d©  = 


\n/2 


=  2^  X?i^Tr(l)  J  Jo(krsin°)sin0  cos 2n+2C  dG 


This  integral  can  he  evaluated  hy  using  Sonine's  first 
formula  (•)  and  one  derives: 


(2.2-31) 


e  r 

— Re  A  * 
2  *  A! 


:  rE  ZI  X  /  \2n+1  0n+l/2  p  . 

sln(krV  'v" 


( ')  Ref  .28  p.373 


-  65  - 


roi^ 


When  kr->*3  (or  r»X),  one  has: 


r  sin  kr  T 

W2(kr) 


\ 


2  sin(kr)cos{kr  -[(1  +  n)/2]n} 


Tt 


Thus ,  when  r  »  A , 


e  r 

(2.2-32)  Re  A1  ^  0 


Hence  this  term  can  be  eliminated  from  Eq.(2.2-23)» 

Let  us  now  consider  the  third  integral  of  Eq.(  2.2-23), 

that  is 


(2.2-33)  LI* 
z 


sincp^  +  E2e3in?2)(E1r  * 


E*  ) 
2r' 


d£ 


This  expression  represents  the  torque  acting  on  the 
antenna  due  to  the  re-radiated  waves.  By  introducing  the  expres¬ 
sions  of  the  field,  by  considering  the  real  part  of  the  integrand 
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o 


and  with  the  help  of  Eq,.(2.2-16)  and  (2.2-17)  one  obtains: 


(2.2-34)  M*  = 

2 


cos^cosoj  sin^sinO  sincpj  ^ 


3in0 


sinfi  d© 


Apart  from  the  sign,  Eq.( 2.2-34)  is  equal  to  Eq.( 1.4-8), 
(the  different  sign  depends  on  the  different  relative  phase  of 
the  currents  on  the  two  half  wave-dipoles  in  the  two  cases) . 

By  using,  for  the  first  term  on  the  right  side  of 

(2.2-23): 


£  r 

(2.2-35)  Re 


A„  = 


e  r 

Re 


(El0sincp1  +  E2€sinc?2)B*.a  flS,' 

£ 


the  same  procedure  as  followed  in  deriving  Eq.( 2.2-27)  from 
(2.2-24)  one  has: 


I 


aTC 


e  r 


(2.2-36)  -^-Re  Ag=  Rei-  *0E0  r2e”i|  cos(^cosO)  slnQ  d©  y 


2n 


X  !  ® 


v0 


ikr  sin©  cosco  .  2 
i  Ysin  co  do: 


sin  cp  dco  > 

J 
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■which  gives: 


e  r 

(2.2-J7)  -|-Re  A2= 


cos(lcr)|  J^Ckr  sinO)cosf^cosO)dO 


Analogously  the  last  term  of  the  sum  on  the  right  side  of  (2.2-23) 


ll  r 


H  r 


(2.2-38)  -^-Re  Ay*  -|-Re  j  (H^cos^costp^  H2cpGOsfi2C°S(?2^-o *- 

4 


becomes; 


r 


(2.2-39)  -2-ReA  »Ee;-t“j2  I  *  r2»*lkr 

2  J  a  O  0 


Tl 


cosf^cosOj  cos  6  dQ 


A  j  e 

'(0 


a2ti 

ikr  sinO  coscp 


coscp  dcp 


which  gives: 
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^or 


(2.2-40)  -J-ReA^ 


“'Vo 


I  E  r 
o  o 


ain(kr) 


I 

h 


■josj^cosG)  J.j (kr  sinO)cost6  d© 


Finally,  substitution  of  Eqs. (2. 2-32)  ,(2.2- 34), (2.2-  JJ)  -and  (2.2-40) 
into  (2.2-38),  gives,  for  r>>X, 


A 


75 


(2.2-41)  =  IIP  +  ^^0I0E0ry  2Cos(kr)j  (kr  sinO)cos^oosO  jd©  - 


h 


Aft 


-  r  sin(kr)|  J^(kr  sinO)cos|^  cosoj  cos  6  dO 

4 


The  first  integral  on  the  right  side  of  Eq.( 2.2-41)  can  be 
performed  by  using  Eq_.(1 .4-10) .  One  has  j 


AH 


(2.2-42) 


J^kr  sin©) cost^c os©J d©  =  - 


2 

£?00s 


11 


71  2  2 

+  k  r 


0 


3* 


’When  kr )}  n  /4 
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(2.2-43)  J.j(kr  sin©) cos cos©j  d©  ~  cos(kr) 


For  the  second  integral  one  can  use  the  formula  (see  Appendix  5); 


(2.2-44) 


r 

1  J^(a  sin©)cos(P  cos©)cos  ©  d©  = 


2 

-  -cos 
a 


a2  2 

P  +a 


i 


2.  cosVI/p^a2, 

2a  ~2 - - 

x  +  P 


2a  Si^Jtr-5? 


„2  2 

P  +a  1 

t  2  ~  JL\V 2 

(jx  +  P 


from  which  one  has: 


-  70  - 


When  kr»w  /4,  Eq.(2.2-45)  becomes: 


(2.2-46) 


J^(kr  sin9)cos 


2 

cos  6  -» 


ain(kr) 


Substitution  of  (2.2-43)  and  (2.2-46)  into  (2.2-41),  gives, 
(for  r  >>X) : 


(2.2-47) 


X 


z 


K1  +  2 

z 


I  E 
o  o 


It  is  to  be  noted  that  the  integrals  of  Eq,.(2.2-41)  can  be  per¬ 
formed  in  an  approximate  way,  by  substituting 


sinO 


for 


ooag  ooaO 
”ilnfi - ' 


This  amounts  to  the  same  as  substituting  the  radiation  pattern 
of  the  elementary  dipole  in  place  of  the  radiation  pattern  of 
the  half  wave  antenna. 
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The  resulting  integrals 


(2.2-48)  E  «*  j  J^kr  sinO)sin  Q  d© 
*0 


2  2 

(2.2-49)  P  =  /  J^kr  sin©) cos  0  sin  Q,  d© 


can  be  performed  by  using  Sonine's  first  formula  (' 


(2.2-50)  E  -Jg  Jy2(kr) 


(2^.51)  t  -  J5/2(to) 


When  r>>X,  one  has: 


(2.2-52) 


Jy2(kr)  -►  cos(kr) 


J5/2(kr)  ♦  -ffip  sin(kr) 


) ,  obtaining 


Substitution  of  (2.2-52)  into  (2.2-50)  and  (2.2-51) 
shows  that  the  approximate  expressions  (2.2-50)  and  (2.2-51) 
of  the  integrals  appearing  in  (2.2-41),  have  values  which,  for 
large  r,  tend  to  the  values  of  the  exact  expressions  (2,2-43) 
and  (2.2-46). 


2.3  -  Evaluation  of  M  . 
“ - - - z 


Let  us  evaluate  by  means  of  Eq^.(  2.2-47). 

The  expression  I  can  be  obtained  by  considering  the 
mean  Poynting  vector  S  of  the  incident  wave,  whose  magnitude  is: 


(2.3.1) 


The  energy  W  absorbed  by  a  single  component  antenna  is 

ct 

given  by 


(2.3-2) 


where  A  is  the  absorption  cross  section  of  a  single  component 
antenna. 

By  considering  the  resistance  in  series  between  the 
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terminals  of  the  antenna,  which  was  assumed  to  he  equal  to  the 
radiation  resistance  R,  one  has 


(2.3-3) 


W 

a 


1 

2 


o 


R 


Prom  (2.3-2)  and  (2.3-3)  one  obtains: 


(2.3-4) 


2  \Te  n 
I  o  o 


I  E 
o  o 


Since  ( ' )  : 


(2.3-5) 


AR  «  Z/k2 


Eq.(2.3-4)  becomes: 


(2.3-6) 


I  0  0 


I  E 


O  0 


4 


W 

a 


U) 


For  the  term  HP  appearing  in  (2.2-47),  one  has,  from  Eq. (2.2-34) 


( ')  see  Ref .8  p.500 
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and  Eq.( 1 .4-1  3) : 


(2.3-7) 


W 

K<  B  _  0.9939 

Z  U) 


where  W  .  is  the  power  re-radiated  by  the  turnstile  antenna.  Since 
the  antenna  is  resonant,  the  scattered  power  is  equal  to  the  ab¬ 
sorbed  power.  Then  =  2W^,  hence 

WWW 

(2. >8)  M  =  4  --  -  1  .9878  --  »  2.0122  -- 

z  u  U)  U) 

Prom  (2.3-2)  one  has: 

(2.3-9) 

and  (2.3-8)  becomes: 

(2.3-10) 


2  W  =  A  S 
a 


H  =  --  1.0061 
z  U) 


S/u  represents  the  flux  per  unit  time  of  the  angular  momentum 
of  the  circularly  polarized  plane  wave  across  the  unit  area  nor- 
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mal  to  the  direction  of  propagation  (')• 

Eq.(2.3-10)  gives  a  simple  relation  between  the  torque 
acting  on  the  antenna  and  the  flux  of  the  angular  momentum  of 
the  incident  waves.  Here  A  is  the  absorption  cross  section  of 
a  single  component  antenna. 

Let  us  now  consider  the  case  of  circularly  polarized 
plane  waves  impinging  on  a  half-wave  dipole  antenna,  that  is 
on  one  only  of  the  two  components  of  the  turnstile  antenna. 

In  this  case  one  obviously  has: 


(2.3-11)  M  =  0 

z 

because  in  this  case  the  antenna  re-radiates  linearly  polarized 
waves . 

Thus  K  is: 
z 


(2.3-12) 


K 

z 


I  E 
0  0 


AS 

u> 


as  can  be  seen  by  putting  E^  =  0  in  the  formulas  of  the  prece¬ 
ding  sections. 

Eq.(2.3-12)  shows  that  the  torque  acting  on  a  single 
half-wave  antenna  is  almost  exactly  the  same  as  the  torque  acting 
on  a  turnstile  antenna.  This  is  due  to  the  fact  that  the  turn¬ 
stile  antenna,  re-radiates  circularly  polarized  waves,  hence 


{ ' )  see  for  instance  Ref. 29 
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angular  momentum 


2.4  -  The  cage  of  the  elementary  dipoles. 


Let  us  now  suppose  that  the  antenna  is  constituted 
by  two  elementary  dipoles  at  right  angles  and  in  quadrature. 

The  incident  field  is  given  by  (2.1-1)  and  (2.1-2)  . 

The  currents  of  the  two  dipoles  are  of  the  type: 

(2.4-1)  I.  =  I  elwt2 

(2.4-2)  Ig  =  i  Iq  e10Jti 

Where  Iq  depends  on  the  incident  field,  and  on  the  absorption  , 
coefficient.  The  dipoles  are  assumed  to  be  tuned  at  the  wave¬ 
length  of  the  incident  radiation. 

The  electromagnetic  field  radiated  by  the  two  dipoles 
is  given  by 


(2.4-3) 


_  V*  i(ut  -  kr)  . 

Eio  '  -  nr  e  ainei 


(2.4-4)  E  -  i  e1'"*  -  kr>  qobO 

2nr 


1 
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(2.4-5) 


tt  I°'>  i(wt  -  kr)  .  A 

%  =  -  Wr  6  sin<?i 


(2.4-6)  Eon  =  i 


I  sZ  . ,  .  .  v 

_  ,  o  i(«t  -  kr) 


2  Q  2hc~ 


sinO, 


(2.4-7) 


I  sZ  . ,  ,  «  * 

_  o  i(ut  -  kr)  . 

E2r  =  fl  e  003e; 

2-nr 


(2.4-8) 


I  s 

,T  •  o  i(wt  -  kr)  .  A 

%  m  1  555  0  alIl02 


where  s  denotes  the  common  length  of  the  dipoles. 

By  introducing  these  expressions  into  (2.2-11)  and  following 
step  by  step  the  procedure  described  in  section  2,  one  obtains 
exactly: 


H  (Is)2 

(2.4-9)  »;  -  -  y  -l— 

and 

/ejT 

(2.4-10)  K  =  M'  +  I-f-2  IsE 
'  z  z  k  o 

Besides,  recalling  that  for  the  dipole  one  has: 
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(2.4-11) 


(2.4.12) 


one  derives  from  (2.3-11) 


(2.4-13) 


which,  introduced  into  (2.4-9),  gives: 


(2.4-14) 


z  U) 


Finally ,  by  expressing  E^  in  terms  of  from  EqL.(2«3-11) » 
using  (2.4-13),  one  derives 


(2.4-15) 


Z  U) 


By  talcing  into  account  (2.3-15),  Eq..(  2.4-15)  becomes 


(2.4-16) 
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In.  this  case  the  coefficient  appearing  in  the  rela¬ 
tion  "between  the  torque  acting  on  the  crossed  dipoles  antenna, 
the  absorption  cross  section  A  and  the  flux  of  the  angular 
momentum  is  exactly  1 .  With  the  aid  of  considerations  similar 
to  those  made  at  the  end  of  the  preceding  section,  it  is  pos¬ 
sible  to  show  that  the  same  torque  acts  on  a  single  elementary 
dipole. 


2.5  -  Numerical  values  and  conclusions . 

Let  a  transmitting  station  be  located  on  the  ground, 
and  radiate  a  power  P,  towards  a  space  vehicle.  The  space 
vehicle  is  assumed  to  carry  a  half-wave  antenna  tuned  at  the 
wavelength  A.  A  torque  results  to  be  applied  to  the  vehicle. 
The  order  of  magnitude  of  the  torque  which  can  be  obtained 
in  this  way,  will  now  be  discussed. 

Let  G-  be  the  gain  of  the  radiating  antenna,  h  the 
distance  from  the  vehicle  to  the  station. 

The  Poynfcing  vector  at  the  vehicle  is  given  by 


(2.5-1) 


S  = 


JPG 

4nh^ 


Let  us  assume  that  the  receiving  antenna  is  placed 
in  a  plane  normal  to  the  direction  of  the  wave.  By  substituting 
(2.5-1)  into  (2.3-12)  and  recalling  that  A  =  0.13  for  the 
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half-wave  antenna,  one  obtains; 


(2.5-2) 


.  M 


0.13 


PGX3 
r.ripiy- 
871  h  c 


For  h  =  1000  km,  X  =  30  m  and  G  =  100  (which,  corresponds  to  a 

parabolic  antenna  of  diameter  D  ~300  m) ,  Eq.(2,.5-2)  gives: 

M  =  1.39  10_17J/W 

In  this  case  a  power  of  lO*1^  W  should  be  radiated,  to  obtain  a 
torque  of  the  order  of  magnitude  of  1  erg.  ' 

If  h  *  100  km,  X  =  30  m  and  G  =  100,  Eq.  ( 2.5-2)  gives 

M  =  1.5  10“15J/W 


and  a  power  of  10  W  should  be  radiated,  to  obtain  a  torque  of 
1  erg.  These  simple  examples  show  that  this  way  of  employing 
circular  polarized  waves  to  control  the  attitude  of  space  vehi¬ 
cles,  is  not  of  pratical  interest.  It  seems  therefore  unneces¬ 
sary  to  examine  different  situations  (for  instance  the  general 
case  in  which  the  vehicle's  antenna  is  not  in  the  plane  of  the 
incident  wave)  or  to  consider  the  effects  of  the  attenuation  by 
the  atmosphere,  of  Faraday  rotation,  or  other  particular  que¬ 
stions. 
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3  -  RADIATING  ANTENNA  UNDER  THE  ACTION  OP  AN  EXTERNAL  FIELD 


3.1  -  Introduction 

A  third  way  of  applying  a  torque  to  a  space  vehicle, 
by  means  of  electromagnetic  radiation  will  he  discussed. 

Let  a  ground  antenna  radiate  linearly  polarized  waves 
towards  a  satellite.  Let  the  satellite  antenna,  on  its  turn  ra¬ 
diate  linearly  polarized  waves.  Denote  by  "the  angle  between 
the  direction  of  polarization  of  the  incident  wave  and  the 
direction  of  the  satellite  antenna. 

The  combination  of  the  fields  radiated  by  the  ground 
antenna  and  by  the  satellite  gives  rise  to  elliptically  pola¬ 
rized  radiation.  The  degree  of  polarization  at  any  given  point 
will  depend  on  the  phase  difference  between  the  two  components, 
on  the  relative  amplitude,  and  on  ^ . 

Let  us  examine  the  case  of  a  thin  half-wave  dipole 
antenna  in  the  free  space,  on  which  the  distribution  of  current 
is  given  as: 


(3-1-1)  I  =  Iq  cos  k^jj  eiw1^ 


where  k^  ■  j-f  and  j is  the  distance  from  the  midpoint. 

A  plane  electromagnetic  wave  is  assumed  to  impinge 
on  the  antenna.  A  reference  system  centered  at  the  center  0 
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of  the  antenna  (Fig. 3.1) ,  is  chosen,  in  such  a  way  that  the 
field  vectors  of  the  plane  wave  are  given  by: 


(3.1-2) 


_  i(wt  -  kz  +  cp)  . 
B  =  e  '  ^  i 


(3.1-3) 


_  o  x(wt  -  kz  +  cp)  . 
3  =  —  e  T  n 

—  c 


Fig.  3.1 
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The  force  P  acting  on  the  element  dp  of  the  antenna 
is  obtained  by  summing  the  forces  acting  on  the  elementary  charge 
p  dj  where  p  denotes  the  linear  charge  density,  and  on  the  current 
element  I  dL, 


(3.1-4) 


T? 


as 


IaI  -i| 


The  torque  acting  on  the  antenna  is 


(3.1-5)  M  =  \  Re  |  Pa n  j  d| 


/-  L 


where  n  is  a  unit  vector,  directed  along  the  antenna  and  L  =  X/4. 
By  introducing  Eq..(  3.1-4)  into  (3.1-5)  one  derives: 


(3.1-6)  M  =  \  Re  h']P  E*  n+(lA#JAn.|d| 


V-L 


3.2  -  Evaluation  of  the  torque  acting  on  the  antenna . 


that : 


The  torque  B,  given  by  Eq_. (3»1-.6)j  is  evaluated,  recalling 


(3.2~1) 


1=  In 
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and  that  by  the  equation  of  continuity  dl/dj  =  -  3p/3t,  One. 
derives,  by  using  Eq_.(  3.1-1): 


(•3.2-2)..  P  -  -  I0  Bin(ic|)e 


ion  t 

1 


Further  one  has: 


(3.2-3)  n  =  cosp  k  +  sin(3  cosa  i  +  sinp  sina  j. 


where  (see  Fig. 3*1)  p  is  the  angle  between  the  antenna  and  the  a 
axis  and  a  the  angle  between  the  projection  of  the  antenna  on  the 
x  y  plane  and  the  x  axis. 

Substitution  of  (3*1-2),  (3.2-*2),  (3*2-3),  (3*1-1), 
(3*2-1)  and  (3*1-3)  into  Eq.(3.1-6)  gives: 
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By  inspection  of  (3*2-4)  it  is  apparent  thatj  if 
4  M  is  an  oscillating  function  v/hose  time  average 
vanishes.  Therefore  we  will  suppose  =  w. 

For  the  integral 


appearing  in  Ecj_.  ( 3. 2-4)  5  one  obtains  after  some  elaboration 


(3.2-5)  S 


cod  p  -  1; 


!%  \ 

:  7i  sin '  -  cospjcosp  +  2cos 


(|cosp 


cos  p  +  1 

cos^p  -  1 


Similarly -.for '•  • 


a  I 


(3*2-7) 


/  cos(k1p' 

J-L  j 


ilc^  cosP 


d£ 


i  k 


iT 

2 


one  obtains 


(3*2-8)  T  =  - 2-- 

cos  8 


7t  sinf2cosp)  +  4  cos  1  71 


-  1 


\  2 


,  cosP 

(  ?QoaP - 1 - 

1  *  1  COS  P  -  1  ! 


By  substituting  (3*2-5)  and  (3*2-7)  into  Eq.( 3*2-4); 


-  86  - 


by  taking  the  real  part,  one  derives 


Ol  B 
o  o 


(3.2-9)  M  =  sincp 


2 w 


sin  p  cosa  eina  Tit 


2  2 

(sin  p  sin  a  -  1)  T  -  cosp  S 


1  + 


+ 


r 


sina  sinp  cosp  T  +  sinp  eina  S 


■where  T  and  S  depend  only  on  a  and  p  and  are  given  by  Eqs.(  3.2-6) 
and  (3.2-8). 

When  p  the  antenna  lies  in  the  plane  of  E  and  B. 

In  this  case  one  has: 

(3.2-10)  3-2  T  -  0 

-and  the  torque  becomes 

cl  E 

(3.2-11)  H  »  -2-2  sincp  sina  k 

w 

When  cp  »  a  ■  w/2  one  has: 
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cl  E 

(3.2-12)  M  =  -_°_2  k 

w 


In  the  case  p  =  0  one  has 

(3.2-13)  S  =  j  T  =  5 

4  4 

and  M  is  given  by 


(3.2-14) 


rccl  E 

H  = - 2 —  sincp  2 

4(u 


3.3-  Numerical  examples  and  conclusions 

If  the  ground  antenna  emits  power  W  ,  with  a  gain 

o 

G,  one  has: 


(3.3-1) 


E 

o 


1  , 
h' 


ZW  G 

g_ 

2n~ 


where  h  denotes  the  distance  between  the  ground  antenna  and 
the  satellite. 

If  W  indicates  the  power  radiated  hy  the  satellite's 

9 

antenna  (which  is  assumed  to  be  tuned  at  the  wavelength  X)  and 
R  the  radiation  resistance,  one  has:  . . .  " 


-  88 


(3.3-2) 


2W 

s 

R" 


I  = 
o 


Let  us  suppose  now  that  the  phase  of  the  satellite  antenna 
'■’is  locked  to 'that  of  the  impinging  wave  with  cp  =  ti/2  and  that  a=|3= 
Thus  Eq.  (  3.2-1  2)  hecomes 


( 3*3-3) 


M  = 


!  ZW  GW 

I  g  S 


If  the  ground  antenna  is  assumed  to  he  parabolic  with 
diameter  D  and  with  uniform  illumination  '(thus  G  =  *V/xs), 
Eq.(3«3-3)  becomes: 


(3.3-4) 


zw  w 

__g_s 

tiR 


Prom  (3*3-4)  one  can  obtain  the  order  of  magnitude  of  the  value 
which  the  torque  can  reach. 

Let  us  assume: 


D  =  100  m,  X  =  30  m,  W  =  100  kW,  W  =  10  W,  Tf:  h  =  300  km: 

g.  s 

M  =  3.4  10“9  J. 
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FOl  ^ 


It  seems  therefore  that  this  third  way  of  applying 
electromagnetic  radiation,  to  control  the  attitude  of  a  space 
vehicle,  cannot  have  practical  importance,  since  the  torque 
is  too  small.  On  the  other  hand,  the  appearance  of  the  power’s 
W  and  W  >  under  the  square  root  sign  leaves  little  hope'-of' 
significantly  improving  the  situation,  even  with  higher  powers 
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CONCLUSION 


In  this  report  the  possibility  has  been  discussed  of 
applying  elliptically  polarized  radiation  to  the  attitude  con¬ 
trol  of  a  space  vehicle . 

In  order  to  get  a  criterion  to  assess  the  effectiveness 
of  the  different  methods  and  procedures,  an  analysis  and  a 
discussion  have  been  made  of  the  different  factors  which  may  di¬ 
sturb  the  orientation  of  the  vehicle  during  actual  flight. 

Typical  orders  of  magnitude  have  been  given  for  the  torques  due 
to  all  these  physical  factors. 

In  the  first  section  an  active  procedure  has  been 
discussed.  The  space  vehicle  sends  out  elliptically  polarized 
radiation  and  is  acted  upon  by  a  reaction  torque.  In  particular, 
the  case  of  a  transmitting  turnstile  antenna  has  been  analysed 
in  detail.  A  simple  relation  was  found  between  the  reaction  tor¬ 
que,  the  radiated  power  and  the  frequency. 

The  power  levels  available  on  space  vehicles  are 
rapidly  increasing  year  after  year.  Reasonable  forecasts  can  be 
made  about  the  power  which  Yvill  become  available  in  the  near  or 
even  far  future.  A  comparison  between  the  torques  obtainable 
with  such  powers  and  the  torques  due  to  external  or  internal 
physical  factors  which  may  disturb  the  orientation  is  not  at  all 
discouraging.  It  seems  therefore  advisable  to  keep  in  mind  the 
possibility  of  this  device  of  attitude  control  in  future  de¬ 
signs,  especially  where  a  continuous  and  small  action  is  to  he 
anticipated  for  very  long  periods  of  time. 
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As  a  by-product  of  this  investigation,  we  point  out 
the  necessity  for  the  future  designer  of  talcing  into  account 
the  torque  produced  by  the  emission  of  elliptically  polarized 
radiation,  which  may  be  undesired  and  represent  a  cause  of 
disturbance  for  the  maintainance  of  a  given  orientation. 

In  the  second  section  the  passive  method  has  been 
investigated,  where  a  ground  station  transmits  circularly  po¬ 
larised  radiation,  which  is  received  by  the  vehicle.  A  relation 
has  been  given  between  the  resulting  torque  and  the  relevant 
parameters,  namely  the  power  of  the  ground  station,  the  gain 
of  the  transmitting  antenna  and  the  absorption  cross  section 
of  the  receiving  antenna.  The  torques  obtainable  turn  out  to 
be  extremely  weak. 

Finally,  in  the  third  section  a  mixed  procedure  has 
been  investigated.  A  ground  station  transmits  linearly  polari¬ 
zed  radiation,  which  at  the  receiving  end  is  utilised  to  control 
the  frequency  and  the  phase  of  the  linearly  polarized  radiation 
emitted  by  the  space  vehicle.  Both  fields,  that  emitted  by  the 
ground  station  and  that  emitted  by  the  space  vehicles  combine 
to  give  an  elliptically  polarized  field.  The  torque  generated 
in  this  way,  by  assuming  reasonable  values  for  the  powers  f 
turns  out  again  to  be  extremely  small,  though  somewhat  better 
than  in  the  previous  case. 

In  conclusion,  it  seems  that  the  active  device  de¬ 
scribed  in  the  first  section  is  both  the  simplest  and  the 
most  effective,  and  the  only  one  to  present  a  practical  inte- 


APPENDIX  1 


In  section  1.2  the  following  relation  has  been  used: 

f  i  r  i| 

(A1-1)  J^div  (u  v)j  r  =  div^(u'r)xjj-  u  v 

where  (u.  v)  is  a  dyadic.  In  order  to  demostrate  Eq_.(A1-1),  let 
us  apply  the  following  formula  ( ' ) : 

(Al-2)  div  (u  v)  =  u  div  v  +  (grad  u) .v 

By  applying  (Al-2)  to  (the  two  parts  of)  (A1-1)>  one- has: 

(u  div  v)  r  +  j~(grad  u)  »v  /  r  =  ,UAr)div  v  +jgrad  (ur.r)J  .v  -  u/v 

By  introducing  rectangular  coordinates  x,  y,  z  (unit  vectors  = 

=  i1 »  i2»  ij,  one  has: 

/  n  v  t)u  .  8u  .  9u  . 

(grad  u)  .v  =  x,  .v  +  3-  i2-I  +  i^v 

thus: 

fcrad  u).vj.-r  =  i2-2  jJ'E*  i3-I  J5'-£ 

( • )  Ref .7  p .26 


-  93  - 


Besides ; 


i-r)  ♦ 


£rad(u-r)  > 


hence: 


(A1-4) 


,  ,~i  3(u/‘r).  3(u^r).  9(u«r). 

grad(u/r)J  .y  =  —  *Z  +  --§“-i2*I  +  — 


fau  3r 

V*  «'■£  +  S'5; 


+  . 


By  making  use  of  (A1-3)  and  (A1-4),  Eq.(AI-l)  becomes: 

.  ,  I  6  \  3r  ( .  \  3r  { .  \  3r  \ 

(A1-5)  0  +(i2-ljuA^  +  (x3.vju/.^  ;  -  fi/v 

but  r  is  the  distance  from  the  origin  of  the  reference  system, 
then: 


Qr  _  .  3r  ,  3r  _ 

"  -1  ’  ffy  “  -2  9  5z  -3 


and  the  term  in  the  parenthesis  in  (A1-5)  becomes: 


v  u/i,  +  v  u/\i„  +  v  u/i,  »  u/.v 
x - 1  y—  —2  z - 3 - 

and  (A1-5) ,  hence  (A1-1)  is  verified. 
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APPENDIX  2 


It  is  to  be  demonstrated  that: 

( div  AU)/1  r  dV  =  0 
'  V 

when  V  is  a  spherical  volume  centered  at  0  (see  p.5  for  the  no¬ 
tations).  The  integral  can  be  written: 


and  because  of  (A1-1): 


where  I  is  the  boundary  surface  of  the  volume  V,  and  n  the  outer 


-  95  - 


normal.  The  last  expression  vanishes  because: 


■  r  (ijj./nJCi^.n)  =  r(n.i1i^n  +  -•~2^2*s+a*-3^3/l-)"0 
k  it 


APPENDIX  3 

In  order  to  perform  the  integration  appearing  in 
Eq.(  1.4-10)  let  us  consider  the  power  series  development  of  J,j(a  sinO) 


co 


> 


(.1) 


m|a  sinO 

\-r" 


i2m+l 


(AJ-1)  J,(a  3in0)  =_n - i:I£-rT7T 


One  has 


(A.$-2)  I  J.jCa  sin'®)  cos(P  cosO)dO 


Am  — ITtmluT —  lsln2m*'e  00a<l3  C03«>ao 


Xs-  („\2o+’ 

*  ^  4-m - TeIt7 —  flj 
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By  putting  m  +  1  *  a,  one  obtains 


(A3-3)  /^(a  sin^)  c°s(P  cosO)dfl 


By  using  Glaisher' s  formula  (')  and  observing  that  =  °* 

one  derives 


U3-4) 


J .  (a  sinO)oos(p  cos£)dD  =  -  -•  cosjyp  +  a 

1  { x 


APPENDIX  4 


Electromagnetic  torque  on  g.  turnstile  antenna,  evaluated  by  means 
of  the  mutual  actions  between  charges  and  currents. 

The  current  distributions  on  the  two  half-wave  antennas 
constituting  the  turnstile  of  Fig.A4-1  are  denoted  by: 


Pig.A4-1 


(«)  see  Ref. 28  p.140  Eq_.(3) 
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(A4-1) 


ii  a  I0  cos(3£y1)  elU)t 

±2  =  i  Io  cos(kz/;))  eiwt  k  =  |- 


Tlie  distributions  of  charges  are: 


o,  ■  r  eiBt 

(A4-2) 

1>2  -  5  Vin(toi2)  eiMt 


Here  and  in  the  sequel  the  convention  for  the  indexes  1  and  2 
is  the  same  as  established  on  section  1.3.  In  (A4-1)  and  (A4-2) 
y1  and  z2  represent  the  distances  of  two  points,  of  the  antenna  1 
and  of  the  antenna  2  respectively,  from  the  origin  0  of  a  system 
of  rectangular  coordinates.  The  two  antennas  cross  each  other  at  0. 
i ,  j_,  k  represent  the  unit  vectors  of  the  reference  system.  The 
time  average  dP  of  the  force  acting  on  the  element  dz^  cf  antenna 
2  at  Q  =(0,0,z„),  due  to  the  distributions  of  charges  and  current 
on  antenna  1  is: 


(A4-3)  dP  =  \  Rej  p2(Q)  E*(Q)  +  i2(Q)k  aB*(Q) 


dz. 


where  1^(Q)  and  B.j(Q)  denote  the  electric  vector  and  the  magnetic 
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induction  vector,  produced  at  Q  by  the  charges  and  currents  on 

antenna  1 . 


Let  and  cp^  represent  the  vector  and  scalar  potentials, 
J.u.0  to  antenna  1.  One  has,  at  a  point  M  =  (x,  y,  a): 


riA/  4 


( A4-4)  A1  =  =  1  j 


.  ko  iwt 


i  (  y  ) 

VM'  -ihr 


4n 


V 


ay- 


■  L  A/4 


(A4-5)  tp1 


o 


fx/4 

iwt  f 


pi(yD 


A/4 


-ikr 

e 


Where:  >  ■  . . . — 

,/  2  ,  s2  2 

r  =  w  x  +  (y  -  y.,)  +  z 

The  field  vectors  at  Q  are: 


(Q)  *  -  gradQ  ^  -  i  w  A^ (Q) 

B1  (Q)  =  outIq  A1 


Since  Q  is  located  on  the  z  axis: 


-  99  - 


(A4-6) 


(A4-7) 


and: 


5z  Hz  4ue 


9  1  iwt  /  -i 


/ 


-  A/4 


..  lf  2  ,  N  2  2 
-lki/s  +(y-y  )  +s 

I^sin(ky t fly- 

2  .  .2  2 
x  +(y  -  yn)  +  z 


c  o  *■  1 


^A/4 


1  iwt  iIo  I  .  /  v  ~ik^:2+ 

i«-  *  -c~| 

0 


,  \  2  2 
-V  +z  h-. 


v  -  A/4 


^7JSg. - aas.  ay 

|  x^+(y-y1 ) ^+s  y 


At  Q,  since  x  and  y  are  equal  to  zero: 


1A/4 


-il 


Si  ^1/Q 


4rce  c 

o 


sin(ky 


,)e-ik^f  ♦  a2  fl  - 


-  A/4 


dy.. 

l  %i~4)  I 


The  function  to  be  integrated  is  odd  with  respect  to  y  ■  0, therefore 


(M-8) 


=  0 
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Thus : 

(u.) 

-  -  H«rJ 

(A4-9) 

and: 

(A4-10) 

/  9?1\ 

I/Q)  =  ~|  5-- 1  -  iw  a^q) 

\  /Q 


Hence  Eq«(A4- 3)  becomes: 


(A4-11)  dP  =  ^Rel 


rA  *  /qa  \* 

"  P2(Qfy~)Q  1  +  iwP2(Q)A1(Q)-i2(Q)^5^  dz2j  = 


“  “  ?Re  P2^%r)  “  i«p2(Q)A*(Q)^i2L^ 


1  dz. 


where  p2,  i 2>  A1  are  evaluated  at  Q* 


One  has: 


\/4 


N,r  io 


( A4-1 2)  -  ?Re  P2J5y-J  =  —--5  ain(kz2)  J.y  oinOiy,) 

0 


sin  (kr) 


ay*. 


»-  X/4 
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,  n  »I 

(A4-13)  ^  Re(ia>p2A*)  - - sin(kz2) 


/A/4 

oos(ky1 )  --~---dy1 

/-X/4 


and; 
(A4-14) 
1 


4  a/ 


**  Re  i  ,  77 — 

2  2  dz 


*  T2 

H  1 

0  o 


/  A/4 


*  “BiT  cos^z2)  35  /  oos(k3ri) 

X/4 


sin(kr)  dy 


Thus  Eq^.(A4-1l)  becomes; 


A/4 


(A4-15)  dP  =  g2  I^Jsin(kz2)  J-  f  sin(ky1  dy1  - 

/-  X/4 


A/4 

-  k  sin(kz2)  |  cos(ky1  dy1  + 

-  A/4 


^A/4 

003(ky1)£i2<S2  jy, 


X/4 


1  dz. 


Let  us  now  consider  an  element  dy^  of  antenna  1 ,  placed  at  the 
point  Pr(0,  y^ ,  0).  The  force  dP*  acting  on  dy.j ,  due  to  the 
charges  and  currents  on  the  antenna  2  is: 
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where  notations  analogous  to  those  of  Eq.(A4-3)  have  been  used: 
One  has: 


i2(2)  =  -  gradp  cp2  -  iwAg(P) 
B2(P)  =  curl  A2 


where: 


A/4 


(A4-17)  A2  *  Agk  =  -2  k  eiu)t 


12<'Z2^  -ikr  , 
-----  e  dz2 


/-  X/4 


( A4— 18)  <p2  -  ^-re 


A/4 

1  iut  f  P2^Z2^  -ikr 

- e  dzn 

0  1  2 
I-  X/4. 


The  following  relations  can  be  derived,  analogous  to 
EqLs.(A4-6)  through  (A4-10): 


(A4-19) 


(»1 


(A4-20) 

/  wTO  - 

( A4-21 ) 

ht\ 

^r/p 

Hence  ; 

( A4-22) 

S2(p)  = 

(A4-23) 

s2(p)  -  -1 

=  0 


/8A2j 

v5?-/P 


I 

;p 


E<q.  (A4-13)  becomes: 

(A4-24)  _ 

i  r  ^ 

«’  -  1  RJ  -P,(p)  BV  s  ♦  ♦  i,c«,  By-,1  >ifay 


A  X 


;  3  A. 


'  S 


Re 


/  Qcp  2- 


L 


pi(p)\^r|p  ■  iuPi(p)vP)  +  vp)i  - dy 


Further  one  has:. 

(M"25)  k  \* 

i  /aV 

“  2  Re  p1(P)'B5"L 


-  I 


/V4 


2-5  sin(ky  )»  I  sin(kz,)-ei5iS£2  te, 
8*e  10 z  ;  2  r 

0  L  X/4 
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( A4-26)  Jse  imp.,  (P)A*(P) 
(A4-27)  -  ^Re 


9  f 

,  V4 

Iun  I 

/ 

X/4 

„  T2 

A  X/4 

^0  0  /, 

f  oos(Rz2)— 5<S£2a22 

V-  X/4 

let  the  points  P  an  Q  he  chosen  at  equal  distances  apart 
from  the  center  of  the  antenna. 

By  comparing  Eq.(A4~25)  with  Eq.(A4-12),  Eq.(A4-26)  with 
Eq .  ( A4- 13),  Eq .  ( A4- 27 )  with  Eq .  ( A4- 14),  and  taking  into  accoun t 
Eqs.(A4-11)  and  (A4-24),  one  has: 

dP.jj,  =  -  dP'  .k 

Pig.A4-2  shows  the  directions  of  the  forces  acting  on  the 


Pig . A4- 2 
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elements  dy.j  and  dz and  indicates  that  the  moments  of  the 
two  forces  with  respect  to  0,  are  of  the  3ame  sign. 

The  integral  in  (A4-15)  will  be  performed  by  means 
of  a  series  approximation  of  the  function  ----- : 

sin  hr  ,  1,3  2  1  ,5  4 

- =  k  -  k  r  +»kr 

r  6  3‘ 

The  function  and  its  derivatives  with  l-espeot  to  y  and  z, 
calculated  at  point  Q,  are  then  approximated  as: 

(«-*»  -  *  - 1  >4? +  h  f  ? ♦  4p5 

(A4-29)  [g-  “  1  k3y1  -  ?T  k5(y2  +  vA) 

(A4-30)  fa  \  ^  k5fj?  ♦  4^2 

By  substituting  (A4-28),  (A4-29)  and  (A4-30)  into  (A4-15) 
and  performing  the  integrations,  one  obtains 

k2!2  £2  ,  2  2 

(A4-31)  dP.i  = - -2&J  z  cos(kz)  +  ^  kz  sin(kz)  - 

4  *eoc  9 

-  z  cos(kz)  -  --^---2  0.635 


. fr 
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Eq_.  ( A4-31 )  gives  the  distribution  of  the  force  acting 
on  antenna  2,  due  to  antenna  1.  Eq.(A4~31)  will  now  he  used  to 
evaluate  the  resulting  moment  T  of  the  distribution  of  force, 
with  respect  to  the  center  of  the  antenna. 

This  moment  is: 


T  = 


A/4 


dP/.Sgk 


A/4 


And  one  obtains  the  result: 

8  2 

(A4-32)  T  =  1  .9  10  XI  i  (II  Joules,  A  meters,  I  Amperes) 


The  torque  H  acting  on  the  turnstile  antenna  is  twice 
the  value  indicated  by  Eq_.(A4-32)j  therefore: 

2 

M  =  3.8  10"8  A  I2  i  =  71.6  L.  i 

—  —  u  — 

By  introducing  the  radiation  resistance  R  =  73*12  ohm,  of  the 
half-wave  antenna,  one  obtains: 

W 

M  =  0.98  -  i 

—  w  - 

where  W  is  the  power  radiated  by  the  turnstile  antenna. 
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This  approximate  result  is  to  he  compared  with  the 
exact  result  given  on  p .  21 . 

M  =  0.99  39  -  i 

03“ 

The  difference  is  practically  negligible. 


APPENDIX  5 


By  following  the  procedure  of  Appendix  39  one  has, 
for  the  second  integral  on  right  side  of  (2.2-41), 


(A5-1) 


I  J^(a  sinO)cos(p  cosQ)cos  0  dO  = 


V, 


0 
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Recalling  Lommel's  formula  (')>  (A5-1)  becomes 


(A5-2) 


J  (a  sinO)cos(P  cos©) cos  0  d©  ® 


(  •)  see  Ref  .28  p.HO  Eq..(l) 


-  109  - 


BIBLIOGRAPHY 


1 

(  )  R.E.  Roberson  'Attitude  Control  of  Satellites  and  Space 
Vehicles'  in:  Advances  in  Space  Sciences,  vol.II,  Aca¬ 
demic  Press,  London  I960. 

2 

(  )  A.  Sadowsky  ’Acta  et  Commendationes  Imp.  Universitatis 
Jurievensis'  7,  N.1-3,  1899,  8  N.1-2,  1900. 

( 'j  P.S  Epstein,  Ann.  d.  Phis.  44,  593,  1914. 

(^)  J.H.  Poynting,  Proc.  Roy.  Soc.  A.  82,  560,  1909. 

(J)  R. A.  Beth,  Phys.  Rev.,  50,  115,  1936. 

6 

(  )  II.  Carrara,  II  Nuovo  Cimento ,  6,  1  ,  1949. 

7 

(  )  G.  Toraldo  di  Francia,  Electromagnetic  Waves , Interscience 
Publishers.  H.Y. 

0 

(  )  R.W.P.  Kin^  ,  The  Theory  of  Linear  Antennas,  Harvard  Uni¬ 
versity  Press.  Cambridge  Mass.  1956. 

(  )  S.A.  Schelkunoff  and  H.T.  Friis,  Antennas  Theory  and 
Practice s John  Wiley  and  Sons,  London  1952. 

("j  R.P.  Ha viland, ' Synchronous  Communications  Satellite'  in: 
Communications  Satellites,  Academic  Press,  London  1962. 

11 

(  )  N.W.  Snyder,  'Power  Supplies  for  Space  Vehicles',  Proc. 

of  XI  I . A. C.  Stock  . olm  i960,  Springer— Ver lag,  Wien  1961. 

1 2 

(  )  Missiles  and  Rockets,  August  15,  24,1960. 

0^)  J .R.  Wetch,  II. M.  Dieckamp  and  G.M.  Anderson,  'The  Practi¬ 
cal  Application  of  Space  Nuclear  Power  in  the  1960's', 
?roc.  of  XI  I. A. C. , Stockholm  I960,  Spring er-Verlag,  Wien 

1961 . 

*1  A 

(  j  Missiles  and  Rockets,  July  23,  12,  1962. 


110  - 


(  '■*)  Missiles  and  Rockets,  September  24,  25,  1962 

(  °)  R.E.  Roberson  'Attitude  Control  of  a  Satellite  Vehicle' 

Proc.VIII  Intern.  Astronaut.  Federation  Congr.,  Barcelona, 
1957,  Springer-Verlag ,  Wien  1958. 

(**')  R.E.  Rober'son  ’Methods  for  the  Control  of  Satellites  and  Spa- 
ce  Vehicles1  vol  I,  Wadd  Techn.  Report  60-643,  July  31,1960. 

1 8 

(  )  R.E.  Fischell,  'Magnetic  and  Gravity  Attitude  Stabiliza¬ 

tion  of  Earth  Satellites'  in  :  Space  Research  II,  pro¬ 
ceeding  of  the  Second  International  Space  Science  Sympo¬ 
sium,  Florence  1961  -  North  Holland  Pub. Comp.  1962. 

(19)  W.R.  Bandecn  and  W.P.  Manger  'Angular  Motion  of  the  Spin 
Axis  of  the  Tiros  I  Meteorological  Satellite  due  to  Ma¬ 
gnetic  and  Gravitational  Torques'  N.A.S.A.  Techn. Note  D— 571, 
April  1961. 

(20)  G.  Colombo  'The  Motion  of  Satellite  1958  Epsilon  around 
its  Center  of  Mass.'S  A  0  Special  Report  No. 70,  July  18, 

1961 . 

21 

(  )  Tatsuso  Obayashi,  'Propagation  of  Solar  Corpuscles  and 

Interplanetary  Magnetic  Fields'  J.  Geophys.  Res.  67,  1717* 

1962. 

(22)  R.H,  Wilson  jr.  'Rotational  Magneto dynamics  and  Steering 
of  Space  Vehicles'  ,  Proc.  XI  I.A.C.,  Stock  iolm  i960, 
Springer-Verlag,  Wien  1961. 

C2-3)  W.G.  Gillespie,  D.G.  Eide '  and  A.B.  Churgin  'Some  Notes  on 
Attitude  Control  of  Earth  Satellite  Vehicles',  N.A.S.A. 

Techn. Note  D-40y  Dec. 1959* 

2A 

(  )  F.L.  Whipple,  'The  Meteoric  Risk  to  Space  Vehicles'  in: 

Vistas  in  Astronautics,  Pergamon  Press,  London  1958. 


Ill 


25 

(  )  W.E.  Frye  and  E.V.B.  Stearns  'Stabilization  and  Attitude 

Control  of  Space  Vehicles'  A.R.S.  J.  No. 12,  927,  1959. 

,26K 

(.  )  Juno  II  Summary  Project  Report, Vol . I, Explorer  VII 

Satellite  N.A.S.A.Techn.Note  D-608,  July  1961 

27 

(  )  H.J.  Pratt  jr.  'Propagation,  Noise  and  General  Systems 

Considerations  in  Earth-Space  Communications'  I.R.E.  Trans, 
on  Comm.  Syst.  vol.CS-8  No. 4  —  221,  Deco1960. 

28 

(  )  G.N.  Watson  'Theory  of  Bessel  Functions'  Cambridge  Univer¬ 

sity  Press,  1958. 

(^)  G.  Toraldo  di  Francia  ,  Boll.  Un.  Mat.  It.  JH,  332  (1956). 


112  - 


VO 

LA 


I  CM 
<4  ft  VO 


CM 

P  LA  H 
o  O  ai 

a  a 

f5  T"  ’H 

-P  VO  Ph 

a  _  _ 


Slg^ 


EH  •  ft 
H  ^  O 

EH  ft  a 
◄  a  R 

Q  ft 

B3-H 

EH  SlO ft 
cj 

ft  o  o 

O  03  ^ 

Hi  S  * 

S".  S 

H  PM  EH 


S  O 
O  H  P 


00  £  0  « 3  P 

ft  Q)  p  d  -rl 

5  Pft  P  3  £ 

§  o  3  a  2  0> 

P  0  0  P  W) 

d  «  ft  u  p  a  d 

4)  0  p  fl>  w  ft 

«>  -h  -p  ;h  s  p 

P  P  0  P  O  d  a 

d  ft  ft  3  ft 

0  0  O  ft  P  ft 

cj  ftp  cj  co  cj 

p  o  P  o  ft  0)  d 

d  ®  ft  >> 

ri 

o  o  p  d  H  £ 

ft  O  ft  EH  ft  a 

ftp  o  d 

cj  p  CQ  *  ft  0)  •* 

ty)  0)  CO  P  0 

ft  >,  J>  0  0  0  H 

P  H  cj  H  d  d  O 

CQ  ft  £  O  S)  O  P 

>  i  -h>  "S 

d  0  0  O  CQ  >> 

p  o  tsj  >  R  U 

p  P  P  P  0) 

H  R  ®  op  o 

a  >>  a  o  o  cj 

O  P  H  Cj  H  o  ft 

p  p  O  ft  o  P  CQ 

P  H  ft  CQ  EH 

CD  p  cfi  Cj 

d  P  >>  ft  H  • 

O  P  H  O  ft  ft  d 

0)  CQ  H  >5  CD  O 

p  CQ  Cj  p  H  CQ 

P  O  O  O  ft  0  ft 

ftP  d  ft  £  0 

'P  p  Cj  o  o 

CD  ft  d  CQ  Cj 

p  P  O  ft  ft  H 

PH  OOP  ft 


CM 

P  LA  H 
OO  n 

g~.S 

P  vo  ft 

£  _ 
O  ft  « 

O  3  EH 


|  CM 

^  ft  vo 

&  •r'l  ® 
H3  6 
Ph  R  ® 
O  o  ft 

CQ 

0 

8  O  Cj 

EH  •  ft 
H  <  O 
EH  •*  d 
EH  ft  a 

£  r 

O  ft 

B  a  .H 

EH  ^0  ft 
Cj 

ft  O  o 
O  cq  ft 

N 

omu 

Ph  <!  O 

en  ft  ^ 

&  * 


u 

o  o 

H  P 

3  8 

S  CO  u 

a  o 

O  H  P 

«  a  o 
p  a  o> 

O  |>  ‘H 

M  o 

ft  ft  R 
ft  O  ft 


H  &S„s 

00  gs  0)  0  ft 
ft  ®  p  £  *h 

Sj'H'jPfi  ^ 
a  o  £  £  ®  ® 

“  £  0  0ft  to 

d  cq  p  r  ft  a  g 

CD  0)  P  CD  ^  w 

o  p  p  p  a  . 

p  P  Cj  p  o  a  *3 
d  p  p  a  *h 
0  0  0  ft  ft  ft 
cj  ftp  0  03  0 

P  o  P  CD  *H  P  R 

d  ®  ft  >» 

a  ft  r  R  o 
o  o  p  R  H  a 

H  0  ft  EH  ft  ® 


ft  0  ft  EH  4  ft  ® 
P  p  OS 

«J  ?  to  «.  -H  ®  * 

t*0  ©  «  +| 

•H  t>j  >  ©  ®  ®  H 

n  a  “  o  |§-H 

f  9«afi ;| 

d  q)  q3  o  p  p 
•h  o  N  >  R  R  « 


p  p 


■S  IIS 


©  O  ‘H  ® 


V/  ^  ^ 

O  -P  H  C3  I-}  ®  £< 

•H  -rl  o  p<  &  £  ® 

■PrlP|8,.Brf 
©  -rl  S?  03 

hp  t>>tH  g  js  H 
o  p  i-h  o  *A  ft  a 

40  ft-H  2  ® 

•iJ  -P  -P  «3  O  O 

<D  ft  d  CQ  Cj 
p  P  O  P  ft  H 
ft  p  o  o  ft  ft 


CM  P 
P  LA  H  d 
o  o  a  a> 
a  d  d 
a  r-  ‘H  ft 
p  VO  Ph  P 

a  _  a, 

o  P4  cj 
o  ^  EH  w 


4i  H  S 
ft  5'. 

P  o 
Ph  R  ® 
OoP 
CQ 

w  a 
§  o  a 

EH  ♦  H 
H  ^  O 
EH  *»  d 

EH  H  a 
<4  d  R 
O  ft 

ft  H  ft 
EH  g)  ft 

O  03  ft 

P  S  ^ 

OpR 

«  •  O 

Eh  ft  Eh 


<H  «  . 

0  I  b  ^ 

H  O  P 

(Do  ^  a)  a  p 

p  0  p  d  *h 

a  p  ft  p  a  * 
do  d  d  2  ® 

p  0  a  p 


d  03  h  r  p  a  a 

0  ©  p  0  a  h 

©  p  p  <h  a  _  p 

p  p  a  p  o  d  cj 

jh  p  p  a  *rl 

CQ  0  0  ft  P  ft 

a  ftp  a  cq  a 

p  OP  0’H  0  R 

R  CD  ft  >> 

d  ft  R  R  Jj  q  ft 

o  o  p  r  h  a 

P  CD  p  EH  ft  a 
p  p  o  d 
a  p  03  c  p  0  *» 
t^o  0)  CQ  P  0 
•ri  >)M®  JH 

ph  a  h  a  a  q 

0  ft  £  O  M  O  P 

©  ft  H  0  p 

>  a  ft  p  d  p  0 
d  0  ^3  0  0  > 
P  O  N  r>  R  R 

p  p  P  H  o 
H  R  0  o  p  o 

a  >>  a  O  0  t  j 

O  P  H  cj  H  o  Pi 
P  P  o  ft  0  &  CQ 

P  H  ft  fH 

0  P  Cj 

r  p  >>p  a  •  . 

o  H  H  o  P  ft  a 
0  0  H  >»  ©  o 
p  0  a  H  H  CQ 

p  o  O  O  ft  0  ft 
ftp  R  ft  a  0’ 
<5  p  p  a  o  o 
0  ft  d  cq  a 

p  p  O  p  *rH  H 

P  H  O  O  ft  ft 


CM  P 

p  la  h  a 

o  o  a  o 


P  VO  fe 

a 


!  CM 
P  VO 

&  p  • 
P  o 
ft  R  ® 
O  o  ft 
0 

0  a 

ft  o  ^ 
p  o  a 

t+  °  *H 
H  <  O 


*<  a  ^ 

O 

B  S  -H 


ft  o  o 
O  0  ft 

o  <i  ^ 

Ofl  h 
«  »  O 

EH  en 


01  o  a 

©  ©  ^  ®  <3  •jj 

rO  O  ?  P  H  'g 

^  O  jP  CQ  0 

s  °5t«  d-p  *> 

d  00  *h  ©  ,a  a  a 

®  ®  P  R  ^  w  T* 

©  p  p  o  a  . 

,0^  ClI'H  5 

d  p  ft  a  p 
0  ©  0  ft  ft  ^ 

a  ftp  a  ®  g 

d  of  ft  0  H 

d  O  ft 

a  ft  r  ?  {3  9  *3 

o  o  r  d  h  8 
ft  ©  p  p  a  ® 

P  p  HOB 


U 

M  O 


a  +3  o  •  d  ®  • 
j  1  | 

d  4  o  o  ^  ^ 

•H  s  .s  ►  ^  ^  © 

H  S»o«2 

«  >»  ®  2  5  «  a. 

O  +=  H  a5  «d  ®  & 

•H-rl  O  P(»^  “ 

p  H  ft  CQ  cH 

Oft  5?  ^ 

d  P  H  JJ  ^ 

q  4(h  pH  O  *h  ft  m 
ffl  03  H  .  >?  S  ° 
cj  ra  d  i-l  H  ®  . 

+3  o  o  o  Qt  ®  ^ 

ft-H  ^  ft?,  ® 
-j;  +»  +»  «J  ®  ® 

§  afj  0  w 

5  “  o  <H  -rj  H 

4»  H  O  O  *d  Pi 


o 

0 

OP 

Cf*  p 

ra  o 
H*  Hj 

cf 

p 

O 

o 

H 

o 

cf 

t?  H 

O 

0 

O 

P 

o 

p 

o 

F* 

CD  O 
0  Oj 

P 

p 

02 

F* 

O  cf 

Fb  0 

o 

o 

H 

O 

cf 

d 

p< 

£ 

cf  <J 

P  0 

F* 

P 

a 

P 

0 

CD 

o 

p 

F* 

<4 

0 

cf 

4 

O 

Q  <+< 
P  P  0 

cf 

p 

cf 

F*  B 

H  >H 

o 

0 

p 

0 

4 

F* 

P 

P 

0 

p 

cf 

Pi 

p 

F* 

B 

cf 

P 

4  F* 

0 

& 

0 

0 

cf 

Pj 

0 

p 

cf  0 

Pi 

4 

0 

P 

pi 

P 

P 

O 

H 

(D  By 

cf  CD 

Pi 

4 

0 

0 

0 

4 

F* 

P 

P 

F- 

Pi 

0 

F* 

(T^ 

P 

0  0 

4 

F‘ 

F* 

Pi 

0 

Pi  P 

O 

0 

B 

4 

0 

0 

CT* 

4 

Pi  0 

O 

0 

ffl 

0 

B 

0 

O 

P 

cf 

P 

O 

O 

F* 

Pj 

cf 

P 

0 

0 

^d 

0 

O 

P 

Fj 

O 

O 

PJ. 

0 

cr  cf 

0 

d* 

Pj 

P 

P 

0 

P 

P 

CD 

cf 

0 

cf 

bi 

d 

d 

0 

0 

F* 

F* 

0 

o 

cf 

d- 

d 

P 

o 

o 

F* 

F* 

02 

B 

B 

O 

0 

cf 

cf 

P 

cf 

F* 

p 

F- 

0 

0 

^d 

H 

B 

O 

0 

cf 

cf 

c/ 

cf 

t.j. 

cf 

P 

P 

d 

p 

0 

P 

O 

0 

O 

P 

CD 

*d 

P 

B 

p 

dj 

P 

p^ 

0 

P 

O 

d 

0 

H 

O 

h 

0 

cf 

<t> 

P 

P 

P 

H 

4 

cf  CD 

H 

o 

H 

O 

<D 

CD 

B 

3 

0 

F* 

P 

O 

F* 

F* 

o  p 

j-j 

0 

P 

F* 

0 

o 

cf 

0 

<1 

o 

0 

4 

»d 

CD 

ESI 

P 

F* 

m 

cf 

CD 

<1 

O 

oJ 

0 

O 

<D 

0 

H> 

P 

^2j 

P 

Oj  02 

B 

O 

O 

CD 

Q 

Fb 

0 

cf  Fb 

4 

P 

PJ 

oJ 

Pi 

H 

cf 

d* 

vD 

d 

rf 

Fb 

4 

P 

P 

cf 

o 

o 

P" 

»Q 

F* 

4 

P 

F* 

p 

W 

cf 

F- 

O 

o 

£•-' 

(jQ 

4 

H)  4 

P 

p 

P 

O 

P 

0 

P 

o 

0 

O 

O 

cf  Fb  4 

0 

P 

cf 

O 

S 

O 

cf 

cf 

0 

P 

H 

H 

0 

cf 

0 

P 

C 

cf 

cf 

CD 

p 

H 

pj 

O 

cf 

B 

0 

0 

0 

0 

F* 

H 

0 

pi 

cf 

cf 

• 

O  cf  P 

0 

0 

0 

0 

P 

B 

P 

0 

O 

P 

P  P 

d-d 

0 

0 

o 

i*  £0 

P 

cf 

8 

cf 

O 

o 

4 

0 

»> 

3  0 

O 

P  0 

P 

cf 

0 

cf 

0 

0 

0 

O 

0 

O 

P 

F* 

P 

cf 

p 

H 

0 

0 

0 

0 

0 

O 

P 

F* 

0 

Fb 

0 

o 

4 

pi 

O 

O 

P 

<4 

4 

0 

O 

4 

Pi 

O 

0 

0 

*2 

P 

cf 

*r* 

CD 

0 

0 

&  0 

Fb 

CD  tO 

P 

• 

O 

F* 

d 

0 

P 

Hi 

H 

P 

P 

P 

p 

*  O 

.r* 

F* 

<1 

0 

p 

0 

O 

P 

FJ 

0 

0 

0 

0 

4 

0 

0 

F* 

O 

O 

P 

4 

0 

0 

0 

Fb 

P 

F- 

0 

P 

e 

0 

C'5 

4 

Pi 

O 

cf 

P 

F- 

0 

K-> 

W 

0 

0 

{.  -! 

4 

P, 

H 

cf  F* 

4 

4 

o 

0 

0 

<4 

ni 

H 

j-j. 

Pj 

H 

cf 

F* 

FJ 

4 

o 

0 

^d 

O 

Q. 

B 

F* 

a 

o 

0 

O 

0 

O 

S 

<5 

o 

p 

B 

H 

g 

B 

cf 

p 

4 

• 

P 

P 

*d 

H 

& 

p 

B 

cf 

4 

o 

cf 

F' 

P 

P 

P 

P 

F* 

F* 

w 

h’ 

cf 

4 

0 

cr 

F* 

P 

0 

Pi 

P 

c 

O 

O 

P 

H 

0 

O 

0 

F9 

cf 

Pi 

P 

6 

O 

O 

p 

P 

0 

o 

F* 

P 

cf 

H 

P 

P 

O 

0 

o 

F* 

P 

cf 

0 

P 

cf 

B 

P 

< 

CD 

O 

0 

4 

cf 

CD 

P 

cf 

0 

0 

p 

cf 

pi 

P 

p 

o 

B 

P 

P 

0 

«- 

0 

>8 

P  cf  pi 

P  O 

P 

P 

F* 

O 

g 

0 

p 

cf 

cf 

0 

P 

P 

0 

P 

F- 

O 

P 

0 

0 

0 

cf 

cf 

4 

4 

4 

o 

p 

0 

0 

o 

F* 

0 

0 

4 

4 

4 

O 

P 

0 

0 

piH 

0 

cf 

O 

P 

P 

p 

O 

P 

a1 

O 

0 

cf 

G 

P 

►d 

P 

P 

o 

O' 

P 

P 

0 

H 

0 

o 

cf 

hd 

P 

p 

O 

P 

P 

0 

O 

P 

p 

0 

P 

P 

0 

P 

p 

p 

0 

O 

p 

P 

0 

0 

0 

4 

0 

P 

H 

cf 

« 

p 

• 

P 

Pi 

4 

4 

0 

0 

H 

cf 

• 

o 

P 

4 

o 

p 

P 

P 

o 

P 

4 

O 

1 

o 

o' 

Pi 

4 

0 

1 

o 

P 

D> 

_ 

_ 

_ 

■  - 

___ 

H 

1 

~ 

0 

_ 

_ 

H 

1 

$ 

H* 

pJ 


O 

P 

0 


H* 

0} 

Cf 

0 

m 

CD 

p 

o 

P 

cf 

P 

0 

a 

F* 

B, 

P 

H 

cf 


CD 

O 

P 

CQ 


CD 

P 

P 


0M  O 
4  CD 


o 

P 


F 
<1 

P  F*  B 
Pj  ^  “ 


m 


CD 


cl-  (D  CO 

0  P 

cf  cl-  p 
F-  CD  CD 

°  P 

P  P 
0 


cf  4 
O 

cf  P 
O  H 

O 

P  4 

CD 

di 


CD 


H 


fit) 


0  d 
H  H  cf 
F*  0  P 
0  O  CD 
cf  CD 
CD  P->  cf 
0  O 

O  4 
cf  P  hQ 
O  p 

$  ® 
4  0  CD 

o  « 

Q 

cf  P  4 
P  F*  F* 
CD  O  CD 
H  P 
<|  Q  O 
CD  **  0 

P  Pi 

F*  3 
O  P  P 
h  a>  ^ 
0  P 
*  P 

P 

4 

0 

-  I  - 


0  _ 
d  h 
*d  p 
H  4 
F*  F* 
0  ESI 
Pi  P 
cf 
cf  F- 
O  O 

P 

cf  • 

P 

0  Jt> 

<4  4 

CD  CD 
P  P 
F‘  O 

O  cf 

H  F- 
CD  O 

•  P 


\3  <+ 
p  o 

0  4 

h 

CD  P 
0  0 
O 

§  $ 
Pi  4 
P 
0 


O  P  CQ  O 
cf  0  F-  Fb 


o 

o  _ 

P  cf  < 
p  P  q 

0  p  d 
CQ  0  4 
Pi  0 
Q 


0 

0  _ 

0  ^  F 
B  a  oo 
BP  P 
Cf  o  H  o 
°  g  3 
B  F-  0 

sf  “  E 

0  dj 
O  O  _ 
cf  Hj  4  p 
O  cf  cf 
O  cf  P  0 
P  B  cf  di 
3  0  a  p  p 
0 


P  cf  0  cf 
-  ■  H. 

cf 

p 

pi 

a> 
o 

cf  Fb 
tf 

P 


& 


cf 

o 

H 

cf 

a 

CM 

O 

0 

o 

o 

p 

P 

4 

0 

F* 

S' 

o 

p 

0 

o 

F- 

P 

d 

0 

0 

P 

< 

0 

P3 

cf 

P 

cf 

P 

F* 

P 

4 

0 

P" 

P 

P 

F* 

P 

0 

F* 

CM 

B 

0 

B 

4 

0 

p 

O 

O 

F* 

P 

cf 

p 

p 

P 

P 

0 

0 

B 

0 

o 

cf 

cf 

cf 

oJ 

cf 

F* 

p 

F* 

0 

p 

0 

P 

O 

0 

O 

p 

CD 

cf 

0 

P 

P 

p 

Si 

F- 

0 

cf 

0 

<1 

o 

0 

4 

O 

CO 

0 

Fb 

P 

d 

4 

p 

P 

P 

P 

H 

►2 

F* 

4 

0 

F* 

0 

P 

cf 

o 

0 

0 

0 

O 

0 

P 

H 

p 

0 

cf 

0 

02 

0 

0 

F* 

P 

0 

Pi 

0 

0 

P 

cf 

0 

cf 

O 

O 

0 

O 

P 

F* 

P 

cf 

P 

O 

4 

P, 

O 

O 

0 

6 

P 

cf 

S 

cf 

0 

p 

0 

P 

0 

P 

0 

0 

Fb 

0 

4 

0 

0 

0 

CM 

4 

Qj 

4 

4 

o 

0 

0 

<1 

di 
p  o 
d  H 
d?  p 
H  4 
F*  F* 
0  W 
Pi  P 
cf 
cf  F- 
O  O 

P 

cf  • 

P 

0  £> 

<1  4 
0  0 


O 

O 

O 

o 

o  p  0  o 

cf 

o 

H 

cf 

F3 

OM 

O 

P 

P 

F* 

0 

Cf  0 

F*  Fb 

0 

o 

O 

P 

p 

4 

0 

0 

cf 

4 

pi  cf  -q 

cf 

F* 

B 

o 

P 

0 

O 

F* 

o 

g 

0  0 

P  cf 

P 

d 

0 

CD 

P 

<i 

B 

cf 

P 

4 

F- 

0  P 

0 

p 

cf 

P 

cf 

d 

H- 

O 

H 

0 

P  d 

Cf  0 

p 

4 

0 

P 

P 

P 

e 

0 

0  4 

F* 

F* 

P 

0 

F 

CM 

0 

P 

4 

P  0 

O  JO 

e 

0 

B 

4 

0 

—  M  *> 

0  O*  cf  0  cf 

0  ^  F-  F* 

B  O  02  cf 

B  0  PP 

cf  0  HO  pi 


O  p 

_  P  F 
P  0  B 
0 


0 


cf  O* 

P  0 

CD  cf 

Si 

cf  0 
O  0 
4  P 
A 

P  cf 
0  P 


0 

_  o 

d>*  cf  H) 

oop 
cf  Hj  4  P  P 

O  cf  cf 

O  cf  P  0  0 

PP  cfd 

^  0  O  PP  d  cf 

0  0  0  O  (D  O 

0  0  4 

F  Ci)  H>  a>  [O 

•  2  &  h1’  m  £ 

o  p  4  o  0  o 

P  F*  0  PJ  p  0 

pj  H  cf  F*  4 

0  O  ^  < 

co  y  f 

P  cf  4  CD 

Pi  P  o 
0  O  F‘ 

Cf  0  P 

Pcf  Pi 
F-O  p  p 
4  4  4  0 


F*  pi 
0 

0  o 
cf  F*  P 
^  “  O 
P 


o 
p 
0 

B 

& 
0 

cf  p  0 

P  P 

cf  cf  P 


P  o 
0 


<4  o 

0  Hj 


cf 
p 

.  0 

H  Pi  0 
F-  Pi 
P 
cf 
P  F* 

pi  O 

P 


V-i 

F*  4 
o  p 


0 


4 
P  di 
Pi  H 
F* 

P 
cf 
0 


Q 


cf 

O 

P 

0 


O 

P 

cf 

„  & 
Fq  o 

Cc  4 
‘  ‘  o 
B 


P 


4 

o 

I 

Pi 


0 

F* 

I 


P 

O  0 
0 

Pi  cf 
O 

o  4 

g  ® 
0  0 


s  g  s 


P*H  0  cf  O 
P  POP 

►d  q  o  p  d  '  q 
4  0  P  H  cf 


cf  cf  0 
0  0  O 

S  3  ^ 

p  p  0 
p  p 


0 
cf  P 
P  F- 

0  O 

^  H 
<4  0 
0  - 
P 


4 

o 


P 

P 

P 


S'  p 

p 

p 

4 

0 


Cf  4 

O 

se 

o-e 

0 

»d 
p  o 
»d  H 
hd  P 
H  4 
F-  F* 
O  ESI 
Pi  P 
cf 
cf  F* 
O  O 

P 

P 

0  > 

<1  4 
0  0 
P  P 
F*  O 
O  cf 
H  F- 
0  O 

*  P 

PJ 

0  4 

0  P 
0  0 
o 

o  cf 

P 

0 


4  0 

0 


